Draft version August 7, 2012 

Preprint typeset using I^T^X style cmulateapj v. 5/2/11 



(N 

o 

(N 
< 

o 

u 

6 



>: 

OV 
O 

in 
in 

o 



X 



A CENSUS OF OXYGEN IN STAR-FORMING GALAXIES: AN EMPIRICAL MODEL LINKING 
METALLICITIES , STAR FORMATION RATES AND OUTFLOWS 

H.J. Zahid 1 , G.I. Dima 1 , L.J. Kewley 1 ' 2 , D.K. Erb 3 & R. Dave 4 

1 University of Hawaii at Manoa, Institute for Astronomy 
2 Australian National University, Research School of Astronomy and Astrophysics 
3 University of Wisconsin Milwaukee, Department of Physics and 
4 University of Arizona, Department of Astronomy 
Draft version August 7, 2012 

ABSTRACT 

In this contribution we present the first census of oxygen in star-forming galaxies in the local 
universe. We examine three samples of galaxies with metallicities and star formation rates at z = 
0.07, 0.8 and 2.26, including the SDSS and DEEP2 surveys. We infer the total mass of oxygen produced 
and mass of oxygen found in the gas-phase from our local SDSS sample. The star formation history is 
determined by requiring that galaxies evolve along the relation between stellar mass and star formation 
rate observed in our three samples. We show that the observed relation between stellar mass and star 
formation rate for our three samples is consistent with other samples in the literature. The mass- 
metallicity relation is well established for our three samples and from this we empirically determine 
the chemical evolution of star-forming galaxies. Thus, we are able to simultaneously constrain the 
star formation rates and metallicities of galaxies over cosmic time allowing us to estimate the mass 
of oxygen locked up in stars. Combining this work with independent measurements reported in the 
literature we conclude that the loss of oxygen from the interstellar medium of local star-forming 
galaxies is likely to be a ubiquitous process with the oxygen mass loss scaling (almost) linearly with 
stellar mass. We estimate the total baryonic mass loss and argue that only a small fraction of the 
baryons inferred from cosmological observations accrete onto galaxies. 



1. INTRODUCTION 

A complete theory of galaxy formation and evolu- 
tion will have to be able to self-consistently account 
for, among other physical processes, the star formation 
and chemical evolution of galaxies. Our understanding 
of galaxy evolution is rooted in the currently accepted 
cosmological model in which large-scale structure in the 
universe traces out the cosmic web of dark matter and 
growth of the universe is accelerated by dark energy. In 
this theoretical framework, a hierarchical formation of 
galaxies is favored in which larger galaxies form as the 
dark matter halos within which they are embedded merge 
over time. It remains uncertain what epoch in cosmic 
history this is the dominant mode of growth. However, 
recent observations of strong correlations observed be- 
tween fundamental galaxy parameters (e.g. mass, age, 
size, luminosity, baryonic content and angular momen- 
tum) have lead some to question the stochastic nature of 
the hierarchical formation scenario (Disney et al. 2008; 
Nair et al. 2010). One possible resolution is that galax- 
ies and groups of galaxies gather matter early on followed 
by quiescent, isolated evolution (Peebles & Nusser 2010). 
The evolution of galaxies may be simpler than a hierar- 
chical formation model suggests. 

A large number of studies have recently revealed that 
there exists a tight relation between stellar mass and star 
formation rates (SFRs) out to z ~ 2 (Noeske et al. 2007b; 
Salim et al. 2007; Daddi et al. 2007; Elbaz et al. 2007; 
Pannella et al. 2009; Elbaz et al. 2011, among others). 
We refer to this as the MS relation. All these studies find 
the slope of the relation to be near unity and a la scatter 
of < 0.3 dex. The relation and its small scatter is taken as 
evidence that secular processes, such as gas accretion, are 



the dominant mechanism for star formation with mergers 
playing a minor role. In particular, Noeske et al. (2007b) 
suggest that the presence of an MS relation with constant 
scatter at several epochs implies that star formation is 
gradually declining with galaxies spending 67% (95%) of 
their star formation lifetime within a factor of ~2 (4) of 
their average SFR. 

Several studies have applied the observational con- 
straints imposed by the MS relation and its evolution 
to uncover star formation histories of galaxies. Noeske 
et al. (2007a) show that their model of "staged" galaxy 
evolution accounts for the observed relation. In their 
model, less massive galaxies have later onset of initial 
star formation with longer timescales of exponential de- 
cay. Similar models result if star-forming galaxies are 
assumed to lie on the MS relation at all epochs. Several 
studies have focused on this simpler approach of conti- 
nuity of star formation along MS relation. Conroy & 
Wechsler (2009) combine this approach with abundance 
matching to dark matter halos, concluding that merg- 
ers play a minor role in mass growth of galaxies. Using 
this approach, Peng et al. (2010) are able to explain the 
shape and evolution of the observed stellar mass function 
for star-forming galaxies. Papovich et al. (2011) apply 
this technique to understand the gas accretion process at 
high redshifts. Leitner & Kravtsov (2011) use this tech- 
nique to show that gas recycling is sufficient to fuel the 
observed star formation in the local universe and Leitner 
(2012) argue that most star-forming galaxies in the local 
universe formed at 1 < z < 2. 

The chemical evolution of the gas-phase of star-forming 
galaxies is largely constrained by observations of the 
mass-metallicity (MZ) relation. Lequeux et al. (1979) 
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were the first to show that the metallicities of galaxies 
increase with stellar mass. The MZ relation is well es- 
tablished in the local universe (Tremonti et al. 2004) and 
has been observed for intermediate (Savaglio et al. 2005; 
Cowie & Barger 2008; Zahid et al. 2011; Moustakas et al. 
2011) and high rcdshift galaxies (Erb et al. 2006a; Man- 
nucci et al. 2009). The shape of the MZ relation is ob- 
served to be relatively constant with evolution in the zero 
point such that galaxies at earlier redshifts are found to 
have lower gas-phase abundance. 

The metal content of galaxies is governed by the pro- 
cesses of star formation and large scale gas flows. Out- 
flowing gas has directly been observed in starburst galax- 
ies (Rupke et al. 2005; Martin 2006; Tremonti et al. 
2007; Rich et al. 2010; Tripp et al. 2011) and is found 
to be a ubiquitous phenomena in higher redshift star- 
forming galaxies (Shapley et al. 2003; Weiner et al. 2009; 
Steidel et al. 2010). A recent survey of the halos of 
galaxies conducted by Tumlinson et al. (2011) reveals 
that large reservoirs of oxygen are found to exist in the 
circum-galactic medium (CGM) of all star-forming galax- 
ies. They conclude that the CGM of star-forming galax- 
ies contains a substantial amount of gas and metals, per- 
haps far exceeding the gas within the galaxies themselves. 
In this study we use our census of oxygen to quantify the 
loss of metals and gas from the ISM of normal local star- 
forming galaxies. 

Census techniques have proven to be crucial in our 
understanding of cosmological evolution. A well con- 
strained inventory of the energy content of the universe 
is one of the greatest triumphs of modern cosmology. By 
comparing a census of the observed baryons in the local 
universe to the expected cosmological density Fukugita 
et al. (1998, among others) showed that the vast ma- 
jority of baryons are not observed. This is one of the 
missing baryons problems. A second related problem is 
that the amount of baryons within galaxies is not in ac- 
cord with expectations inferred from the properties of 
the dark matter halos within which they are embedded 
(e.g Bell et al. 2003b). Theoretical cosmological mod- 
els suggest that the missing baryons are to be found in 
the warm-hot intergalactic medium (WHIM) (Cen & Os- 
triker 1999; Dave et al. 2001). Baryons in this phase may 
or may not be associated with galaxies and it remains un- 
clear what fraction of the baryons accreted onto galaxies 
and were later ejected. 

In this study we present a self-consistent, empiri- 
cally constrained census model of oxygen in star-forming 
galaxies. The data used in this study are presented in 
Section 2 and the methods used in deriving the stellar 
masses, metallicities and SFRs of galaxies are discussed 
in Section 3. We parameterize the SFRs of galaxies as a 
function of stellar mass and redshift in Section 4 by ex- 
amining the MS relation at several redshifts. We describe 
the various components of our oxygen census in Section 5. 
We develop our self-consistent empirical models for the 
star formation and chemical history of galaxies in Section 
6 and 7, respectively, by imposing the continuity condi- 
tion that galaxies build up their stellar mass by evolv- 
ing along the empirical relation between stellar mass and 
SFR with the metallicity inferred from the MZ relation 
at several redshifts. In Section 8 we present the results of 
our census and in Section 9 we discuss systematics and 
uncertainties in our approach. We provide a discussion 



of our results in Section 10 and a summary of our results 
in Section 11. For this study we adopt the standard cos- 
mology (H ,n mi n A ) = (70 kms- 1 Mpc' 1 , 0.3, 0.7). 

2. THE DATA AND SAMPLE SELECTION 

In this section we describe the local sample of galax- 
ies from SDSS (Section 3.1), our intermediate redshift 
sample from DEEP2 (Section 3.2) and a high redshift 
sample from Erb et al. (2006a, Section 3.3). Galaxies 
are primarily selected such that the chemical properties 
can be determined from their spectra. The binned data 
are given in a Table 1. 

2.1. The SDSS Sample 

We draw our local sample from the SDSS DR7 which 
consists of ^900,000 galaxies spanning a redshift range 
of < z < 0.7 (Abazajian et al. 2009). The survey 
has a Petrosian limiting magnitde of rp = 17.8 covering 
8,200 deg 2 . The spectra have a nominal spectral range 
of 3900 - 9100A and a spectral resolution of R - 2000. 
We make use of the ugriz-b&nd photometry available 
for each object (Stoughton et al. 2002) and the publicly 
released emission line fluxes measured by the MPA-JHU 
group 1 . We refer to the sample presented here as the 
"metallicity selected SDSS sample" . 

We correct for dust extinction in the emission lines by 
inferring a reddening correction from the Balmer decre- 
ment. For case B recombination with electron temper- 
ature T e = 10 4 K and electron density n e = 10 2 cm~ 3 , 
the intrinsic Ha/H/3 ratio is expected to be 2.86 (Ostcr- 
brock 1989). We get the intrinsic color excess, E(B— V), 
and the correction for dust attenuation using the extinc- 
tion law of Cardelli et al. (1989) and a corresponding 
R v = 3.1. We note that the results of this study are not 
dependent on our choice of a particular extinction law. 

We select a pure star-forming sample of local emission 
line galaxies from the SDSS DR7. We first distinguish 
star-forming galaxies from AGN by constraining the ion- 
izing radiation source using the [OIH1A5007, [NII]A6584, 
H/? and Ha emission lines (Baldwin et al. 1981; Veilleux 
& Osterbrock 1987; Kauffmann ct al. 2003; Kewley et al. 
2006). In particular, we remove galaxies using the equa- 
tion given in Kewley et al. (2006) where 

log([OIII]/H/3) > 0.61/ (log([NII]/Ha) - 0.05) + 1.3. 

(1) 

In order to avoid aperture effects, we require a g-band 
fiber aperture covering fraction > 30% in addition to im- 
posing a lower redshift limit of 0.04 (Kewley et al. 2004). 
The median covering fraction for the metallicity selected 
SDSS sample is 38%. Kewley et al. (2006) find that the 
SDSS sample is incomplete at higher redshifts and in or- 
der to minimize evolutionary effects we also impose an 
upper limit redshift cutoff of z = 0.1. 

In order to establish comparable samples, galaxies 
in the local sample are selected from SDSS using the 
same selection criteria as the DEEP2 sample. In par- 
ticular, galaxies are selected to have S/N of H/3 > 3, 
c_R23 < 2 and equivalent width of H/3 > 4A. Here, 
&R23 is the error in the R23 parameter which is the ratio 
of the oxygen nebular emission ([OIL1A3727 doublet and 

1 http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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[OIII]A4959, 5007) to H/3. These particular selection cri- 
teria gives us a sample of ^20,000 star-forming galaxies 
in the limited redshift range of 0.04 < z < 0.1. 

2.2. The DEEP2 Sample 

Our sample of intermediate redshift star-forming 
galaxies is taken from the DEEP2 survey (Davis et al. 
2003). Details of sample selection and properties are 
given in Zahid et al. (2011), here we summarize the data 
selection. The survey consists of ^45,000 galaxies tar- 
geted mostly in the redshift range of 0.7 > z > 1.4 by 
applying a color preselection using _Bi?i-band photom- 
etry (Coil et al. 2004). The survey has a limiting mag- 
nitude of Rab = 24.1 and covers 3.5 deg 2 . The spectra 
have a nominal spectral range of 6500 - 9100A and a reso- 
lution of R <~ 5000. The emission line equivalent widths 
are measured in Zahid et al. (2011) and we adopt the 
same values here. 

The sample selection is based on the spectral and pho- 
tometric properties of the galaxies. We reduce AGN con- 
tamination by first requiring that log(i?23) < 1 (Kcw- 
ley et al., in prep). Using the color separation for blue 
and red galaxies parameterized by Willmer et al. (2006), 
Weiner et al. (2007) conclude that only a small fraction of 
blue galaxies in DEEP2 appear to harbor AGN whereas a 
large fraction of red emission line galaxies show evidence 
of AGN emission. We further limit AGN contamination 
by removing 48 galaxies in the sample that are classified 
as red galaxies using the color division of Willmer et al. 
(2006). 

Given the nominal spectral coverage, the redshift range 
of galaxies in our sample is limited by the necessity 
to simultaneously observe both the [OII]A3727 doublet 
and the [OIII]A5007 emission lines which are required 
for chemical analysis. We further require that the S/N 
H/3 > 3, the error of the i?23 emission line ratio, o\r 2 3, 
be less than 2 and that equivalent width of H/3 > 4A. Fi- 
nally due to ambiguity in the metallicity determination 
of DEEP2 galaxies at low stellar masses (see Section 3.2) 
we also remove all galaxies with M* < 10 9 2 M Q . This 
selection criteria gives us a sample of 1348 star-forming 
galaxies in the limited redshift range of 0.75 < z < 0.82. 

2.3. The E06 Sample 

Erb et al. (2006a, E06 hereafter) determine the MZ re- 
lation from 87 star- forming galaxies at z <~ 2.2 selected 
from a larger sample of 114 galaxies described in Erb 
et al. (2006b). The galaxies are selected on the basis of 
their rest-frame UV colors and redshifts are determined 
from their UV spectra. All galaxies in the sample have 
UGRK-hanA photometry. Most also have J-band pho- 
tometry and 32 galaxies have been observed at 3.6, 4.5, 
5.4 and 8.0 /im with IRAC onboard the Spitzer Space 
Telescope. Ha spectra were obtained using NIRSPEC 
on the Keck II telescope. 

The metallicity for these galaxies is determined from 
the emission line ratio of [NII]A6584 to Ha. The 
[NII]A6584 line is sensitive to metallicity with the 
strength of the line decreasing with decreasing metal- 
licity. The S/N of the individual galaxy spectra is in- 
sufficient to measure the weak [NII]A6584 line. In order 
to increase the S/N of their spectra and to increase the 
chance of detecting [Nil] emission line at low metallicites, 



E06 stack 14 or 15 individual galaxy spectra binned by 
stellar mass into 6 composite spectra. The [NII]A6584 
and Ha emission line fluxes and metallicities are mea- 
sured from these composite spectra. 

3. METHODS 

In this section we discuss our methods for determining 
stellar masses (Section 4.1), metallicities (Section 4.2) 
and SFRs (Section 4.3). 

3.1. Stellar Mass 

We use the Le Phare 2 code developed by Arnouts, S. & 
Ilbert, O. to estimate the galactic stellar mass. This code 
estimates the stellar masses of galaxies by comparing 
photometry with stellar population synthesis models in 
order to determine the mass-to-light ratio which is then 
used to scale the observed luminosity (Bell et al. 2003b; 
Fontana et al. 2004) . We synthesize magnitudes from the 
stellar templates of Bruzual & Chariot (2003) and use a 
Chabrier (2003) IMF. The 27 models have two metallic- 
ities and seven exponentially decreasing star formation 
models (SFR cx e -*/ T ) with t = 0.1,0.3,1,2,3,5,10,15 
and 30 Gyrs. We apply the extinction law of Calzetti 
et al. (2000) allowing E(B-V) to range from to 0.6 and 
the stellar population ages range from to 13 Gyrs. The 
median statistical error for the derived stellar masses, de- 
termined from propagating the uncertainty in the pho- 
tometry, is 0.15 dex. Though systematic effects may be 
large (Drory et al. 2004; Conroy et al. 2009), we have 
consistently measured the stellar masses for our differ- 
ent samples giving us a robust relative measure (Swindle 
et al. 2011, R. Swindle, private comm.). 

We account for emission line contributions by tak- 
ing the Kennicutt (1998a) relation between the synthe- 
sized UV luminosity and SFR and the emission lines. 
This treatment accounts for Ha, H/3 and [OIIJA3727 and 
[OIII]A4959,5007 (Ilbert et al. 2009). In our sample, 
making no correction for the emission line contributions 
does not significantly alter our mass determinations. We 
adopt the median of the mass distribution and take the 
68% confidence interval as a measure of the error. In 
Zahid et al. (2011) we compare this method with the 
method used by the MPA/JHU group to determine stel- 
lar masses of the SDSS galaxies. We find that the two 
estimates differ by a constant offset of ^0.2 dex and that 
the dispersion between the two methods is 0.14 dex. 

E06 measure stellar masses using a similar method of 
comparing photometry with stellar population synthesis 
models (Erb et al. 2006c). However, an important differ- 
ence is that they measure the "total" stellar mass, which 
is the integral of the SFR over the lifetime of the galaxy. 
They find that this stellar mass is ~ 10 — 40% higher than 
the instantaneous (or what they term "current living") 
stellar mass. However, most mass estimates found in the 
literature are generally the instantaneous and not "to- 
tal" stellar mass. We recalculate the stellar masses for 
their individual galaxies and bin the data according to 
their original binning. In Figure 1 we compare the mass 
estimates for each of their 6 bins. The x-axis shows our 
mass estimate and y-axis is the original mass estimate 
from E06. The error bar represents the rms dispersion of 

2 http : / /www.cf ht.hawaii.edu/~arnouts/LEPH ARE/cfht_ 
lephare/ lephare.html 
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TABLE 1 
Data 



log(M*/M ) 


12 + log(0/H) 


E(B— 


V) 


log(tf) 


SDSS 


8.51 


8.707 ± 0.004 


0.11 


± 


0.02 


-0.20 


± 


0.02 


8.82 


8.736 


± 


0.006 


0.13 


± 


0.02 


-0.12 


± 


0.02 


8.97 


8.787 


± 


0.007 


0.15 


± 


0.02 


-0.19 


± 


0.02 


9.08 


8.819 


± 


0.008 


0.16 


± 


0.02 


-0.20 


± 


0.02 


9.17 


8.859 


± 


0.007 


0.18 


± 


0.02 


-0.19 


± 


0.02 


9.23 


8.875 


± 


0.006 


0.19 


± 


0.01 


-0.14 


± 


0.01 


9.30 


8.900 


± 


0.006 


0.20 


± 


0.01 


-0.12 


± 


0.01 


9.36 


8.920 


± 


0.006 


0.22 


± 


0.01 


-0.10 


± 


0.01 


9.41 


8.923 


± 


0.006 


0.23 


± 


0.01 


-0.07 


± 


0.01 


9.45 


8.946 


± 


0.006 


0.24 


± 


0.01 


-0.06 


± 


0.01 


9.49 


8.947 


± 


0.006 


0.24 


± 


0.01 


-0.03 


± 


0.01 


9.54 


8.969 


± 


0.006 


0.27 


± 


0.01 


-0.02 


± 


0.01 


9.57 


8.977 


± 


0.004 


0.26 


± 


0.01 


-0.01 


± 


0.01 


9.61 


8.993 


± 


0.006 


0.27 


± 


0.01 


0.03 


± 


0.01 


9.64 


8.989 


± 


0.005 


0.28 


± 


0.01 


0.05 


± 


0.01 


9.68 


9.007 


± 


0.004 


0.28 


± 


0.01 


0.07 


± 


0.01 


9.71 


9.010 


± 


0.004 


0.30 


± 


0.01 


0.09 


± 


0.01 


9.75 


9.022 


± 


0.004 


0.30 


± 


0.01 


0.11 


± 


0.01 


9.78 


9.035 


± 


0.003 


0.32 


± 


0.01 


0.13 


± 


0.01 


9.81 


9.037 


± 


0.004 


0.33 


± 


0.01 


0.16 


± 


0.01 


9.85 


9.048 


± 


0.003 


0.33 


± 


0.01 


0.17 


± 


0.01 


9.88 


9.056 


± 


0.003 


0.36 


± 


0.01 


0.22 


± 


0.01 


9.92 


9.059 


± 


0.003 


0.36 


± 


0.01 


0.24 


± 


0.01 


9.95 


9.068 


± 


0.003 


0.38 


± 


0.01 


0.27 


± 


0.01 


10.00 


9.061 


± 


0.003 


0.39 


± 


0.01 


0.32 


± 


0.01 


10.04 


9.081 


± 


0.002 


0.40 


± 


0.01 


0.34 


± 


0.01 


10.09 


9.084 


± 


0.003 


0.43 


± 


0.01 


0.40 


± 


0.01 


10.15 


9.088 


± 


0.002 


0.44 


± 


0.01 


0.47 


± 


0.01 


10.24 


9.086 


± 


0.003 


0.47 


± 


0.01 


0.53 


± 


0.01 


10.39 


9.095 


± 


0.002 


0.52 


± 


0.01 


0.73 


± 


0.01 


DEEP2 


9.25 


8.69 


± 


0.02 


0.17 


± 


0.03 


0.35 


± 


0.03 


9.32 


8.76 


± 


0.02 


0.19 


± 


0.02 


0.32 


± 


0.02 


9.39 


8.78 


± 


0.02 


0.20 


± 


0.03 


0.37 


± 


0.03 


9.44 


8.77 


± 


0.02 


0.21 


± 


0.03 


0.39 


± 


0.03 


9.49 


8.74 


± 


0.01 


0.21 


± 


0.03 


0.43 


± 


0.03 


9.56 


8.80 


± 


0.02 


0.23 


± 


0.03 


0.49 


± 


0.03 


9.64 


8.83 


± 


0.02 


0.24 


± 


0.04 


0.54 


± 


0.04 


9.72 


8.84 


± 


0.02 


0.26 


± 


0.03 


0.60 


± 


0.03 


9.79 


8.86 


± 


0.02 


0.27 


± 


0.05 


0.60 


± 


0.05 


9.87 


8.92 


± 


0.02 


0.31 


± 


0.05 


0.72 


± 


0.05 


9.97 


8.94 


± 


0.02 


0.33 


± 


0.03 


0.77 


± 


0.03 


10.07 


8.93 


± 


0.02 


0.35 


± 


0.03 


0.82 


± 


0.03 


10.18 


8.96 


± 


0.01 


0.38 


± 


0.03 


0.92 


± 


0.03 


10.33 


9.00 


± 


0.01 


0.43 


± 


0.02 


0.98 


± 


0.02 


10.59 


9.04 


± 


0.01 


0.53 


± 


0.04 


1.29 


± 


0.04 


E06 


9.14 


<8.55 


0.14 


± 


0.09 


1.40 


± 


0.14 


9.56 


8.72 


± 


0.07 


0.21 


± 


0.07 


1.41 


± 


0.05 


9.89 


8.82 


± 


0.06 


0.28 


± 


0.08 


1.60 


± 


0.06 


10.12 


8.86 


± 


0.06 


0.33 


± 


0.07 


1.52 


± 


0.08 


10.32 


8.92 


± 


0.06 


0.39 


± 


0.09 


1.78 


± 


0.06 


10.73 


8.97 


± 


0.05 


0.52 


± 


0.06 


1.96 


± 


0.06 



Note. — The stellar mass (column 1), mctallicity determined 
from the diagnostic of Kobulnicky & Kewley (2004, column 2), fit- 
ted E(B— V) from Equation 6 (column 3) and SFR(column 4) for the 
SDSS, DEEP2 and E06 samples. For the E06 sample the metallicity 
has been converted from Pettini & Pagcl (2004) to Kobulnicky & Kew- 
ley (2004) using the conversion constants in Kewley & Ellison (2008). 
DEEP2 and E06 sample we have corrected for dust extinction when de- 
termining SFRs from the Balmcr lines using the E(B— V) values given 
in column 3. For the DEEP2 and SDSS data the values are the median 
in bins of stellar mass. The errors are determined from bootstrapping 
and are analogous to the standard error on the mean. For the E06 
sample, the errors are the standard error of the mean. Each data bin is 
equally populated such that the SDSS, DEEP2 and E06 samples con- 
tain ~ 700, ~ 90 and ~ 15 galaxies in each bin, respectively. Electronic 
version available upon request. 
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Fig. 1. — The stellar mass determined by Erb et al. (2006c) plot- 
ted against our determination using Le Phare. The dashed line is 
the one-to-one agreement and the solid line is offset by 0.29 dcx. 
In the five higher mass bins, the stellar mass estimates used in E06 
are greater by a factor of 2 (0.29 dcx). 

stellar masses in each mass bin. For the five higher mass 
bins, we find that there is a constant offset such that the 
E06 estimates are 0.29 dex higher than our estimate. For 
the lowest mass bin, there appears to be a near one-to- 
one agreement, though the 0.29 dex offset is within the 
errors. This is most likely due to the fact that the low- 
est mass galaxies are younger and therefore have a much 
closer agreement between their instantaneous and "to- 
tal" stellar mass. Moreover, we do not concern ourselves 
too much with the lowest mass bin as the metallicity 
measure is only an upper limit for this bin. Since we do 
not rebin the E06 data in determining the metallicity we 
adopt their stellar mass values but subtract 0.29 dex to 
make them consistent with our estimates. 

Figure 2a shows the distribution of the stellar masses 
for all three samples. 

3.2. Metallicity 

We use two strong-line methods to determine the 
metallicity of galaxies in our sample. In this section we 
only present the parameterization of the calibration and 
defer a detailed discussion of the methods until Section 
8.3. For the SDSS and DEEP2 sample we determine 
metallicities using the calibration of Kobulnicky & Kew- 
ley (2004). This method relies on the R23 and 032 pa- 
rameters which arc defined as 



i?23 



[OII]A3727 + [OIII]A4959, 5007 

H/3 



and 



(932 = 



[OIII]A4959,5007 
[OII]A3727 



(2) 



(3) 



Here the ratio implies the ratio of the measured line in- 
tensities. We have used the assumption that the ratio of 
the fluxes of [OIII]A5007 to [OIII]A4959 is 3 (Osterbrock 
1989). Due to the higher S/N of the [OIII]A5007 line, 
for the sum of the [OIII]A4959 and [OIII]A5007 flux we 
adopt a value of 1.33 times the [OIII]A5007 flux. 

The R23 strong-line method is known to be sensitive 
to the ionization parameter. The ionization parameter 
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Fig. 2.— A histogram of the a) stellar mass, b) metallicity, c) SFR (M Q yr _1 ) and d) the fitted E(B-V) for the DEEP2 (solid black) 
and SDSS (dashed blue) samples. The values for the 6 binned data points of E06 are shown by the red ticks. 



characterizes the ionization state of the gas and quan- 
titatively represents the number of ionizing photons per 
second per unit area divided by the hydrogen density. 
The ionization parameter can be constrained by measur- 
ing the ratio of the line intensity of the same element 
at two ionization stages. The method of Kobulnicky & 
Kewley (2004) does this explicitly by using the 032 line 
ratio. Because both the metallicity and ionization pa- 
rameter are interdependent, an iterative scheme is used, 
the details of which are provided in the appendix of Kew- 
ley & Ellison (2008). In particular we use equations A4 
and A6. 

The DEEP2 data are not flux calibrated so we use 
the line equivalent widths with a correction applied for 
Balmer absorption (Zahid et al. 2011). For consistency 
we determine the SDSS metallicities using line equiv- 
alent widths as well. Several studies have established 
that line equivalent widths can be substituted for line 
fluxes when measuring metallicities if data are not flux 
calibrated or a reliable reddening estimate is unavailable 
(Kobulnicky & Phillips 2003; Moustakas et al. 2010; Za- 
hid et al. 2011). In particular, we test this on our SDSS 
sample and find that the dispersion between the metal- 
licities measured using equivalent widths and dereddened 
line fluxes is ^0.05 dex which is less than intrinsic un- 
certainties of the strong-line method. The average differ- 
ence in the MZ relation derived using the two methods 
is < 0.01 dex. 

One issue with using the i?23 method is that metallic- 
ity is not a monotonic function of i?23. For a particular 
value of the i?23 parameter there are two possible val- 
ues of metallicity known as the upper and lower branch 
metallicity. This degeneracy is generally alleviated by 
employing a second metallicity diagnostic that relies on 
line ratios that are monotonic with metallicity. For our 
SDSS sample of galaxies, the use of [Nil] /Ha reveals that 
the vast majority (~99%) of galaxies are on the upper 
branch. For our DEEP2 sample, the nominal wavelength 
coverage of the spectra does not allow us to simultane- 
ously observe the i?23 lines and [Nil] /Ha. We assume 
all galaxies with M» > 10 9,2 Mq lie on the upper metal- 
licity branch (for more details see Figure 5 in Zahid et al. 
2011). 

The calibration of Kobulnicky & Kewley (2004) is 
based on theoretical photoionization models. Empirical 
methods rely on calibrating strong-line ratios to metal- 
licities determined using temperature sensitive auroral 



lines. E06 determine metallicities using the empirical 
calibration of Pettini & Pagel (2004). The metallicity is 
given by 



12 + log(0/H) = 8.90 + 0.57 x N2 



(4) 



where N2 = log([NII]A6584/Ha). Using - 28, 000 galax- 
ies from the SDSS DR4, Kewley & Ellison (2008) derive 
constants for converting between various diagnostics. We 
use these conversion constants to consistently determine 
metallicities for all three samples using both the theoret- 
ical and empirical calibrations of Kobulnicky & Kewley 
(2004) and Pettini & Pagel (2004). 

The metallicity of a galaxy is traditionally defined as 
a number density of oxygen relative to hydrogen and is 
given as 12 + log(0/H). For this study we require the 
mass density of oxygen relative to hydrogen. We convert 
number density to mass density using the relation 



\og(Z g ) = 12 + log(0/H) - 



(5) 



12 - log 



/ Mq/M h 
\XM h + YM H( 



:log(0/H)-log 



15.999/1.0079 



0.75 -1.0079 + 0.25 -4.0026 



For the remainder of the paper, Z g refers to the gas-phase 
mass density of oxygen relative hydrogen. Also, when 
referring to metallicity or gas-phase oxygen abundance 
we mean to refer to Z g . 

Figure 2b shows the distribution of the metallicities for 
all three samples. 

3.3. Star Formation Rate 

The SFRs for the SDSS sample are derived by the 
MPA/JHU group using the technique of Brinchmann 
et al. (2004) with additional improvements given by 
Salim et al. (2007). The strong emission lines of each 
galaxy are fit using the nebular emission models of Char- 
lot & Longhetti (2001). All emission lines contribute 
to constraining the dust attenuation, though the largest 
contribution comes from the Ha/H/3 ratio. At the me- 
dian redshift of the SDSS, about a 1/3 of the galaxy light 
is contained within the 3" fiber. In order to account for 
losses, Brinchmann et al. (2004) derive an aperture cor- 
rection which was improved upon by Salim et al. (2007). 
We convert from a Kroupa to Chabrier IMF by dividing 
by 1.06. 



6 




8-0 8.5 9.0 9.5 10.0 10.5 11.0 -0.50 -0.25 0.00 0.25 0.50 

log(tf,/« ) AE(B-V) 



Fig. 3. — A fit to E(B— V) as function of stellar mass and metal- 
licity. The left panel shows the E(B— V) sorted into 8 color-coded 
bins of metallicity and plotted as function of stellar mass. The 
underlying curves are the fitted relation given by Equation 6 and 
color-coded by metallicity. The right panel shows a histogram of 
the residuals between the observed E(B— V) and those determined 
from the fitted relation. 

We estimate the SFRs for the DEEP2 and E06 sam- 
ple from the H/3 and Ha emission lines, respectively. In 
order to reliably estimate the SFR, it is necessary to 
make a correction to the Balmer line flux due to dust 
extinction. Since the nominal wavelength coverage of 
both the DEEP2 and E06 spectra do not allow us to 
observe multiple order Balmer lines, we are not able to 
derive a correction from the Balmer decrement. Garn & 
Best (2010) have shown that for a sample of star-forming 
galaxies from the SDSS the extinction determined from 
the Ha/H/3 ratio can be predicted from the stellar mass, 
Ha luminosity or metallicity of the galaxy with a disper- 
sion of ~0.1 dex. The relation between stellar mass and 
dust extinction found by Garn & Best (2010) is shown to 
be valid out to z ~ 1.5 (Sobral et al. 2011) and a relation 
between dust extinction, stellar mass and metallicity is 
also observed in galaxies at z ~ 2 (Reddy et al. 2010). 

Xiao et al. (2012) obtain a slightly better fit (dispersion 
of 0.07 dex) by incorporating galaxy inclination and us- 
ing a different parameterization. We determine E(B— V) 
for our sample of SDSS galaxies from the Balmer decre- 
ment assuming the Cardelli et al. (1989) extinction law 
and parameterize E(B— V) as a function of stellar mass 
and metallicity using a similar formulation to Xiao et al. 
(2012). The color excess is given by 

E(B-V) = (po+PiZP>) x M ps , (6) 

where Z = io( 12 + lo g(°/ H )- 8 ) and M = M*/10 10 . The 
metallicities are derived using the calibration of Kobul- 
nicky & Kewley (2004). We perform a non-linear least 
square fit using the MPFIT 3 set of routines (Markwardt 
2009). The data are weighted by their la statistical un- 
certainties and the errors are propagated to the parame- 
ters. The fit returns p = 0.12 ±0.01, pi = 0.041 ± 0.006, 
p 2 = 0.77±0.06 andp 3 = 0.240± 0.002. The left panel of 
Figure 3 shows the mean E(B— V) binned by stellar mass 
and color coded by metallicity. The fit to these binned 
data are shown by the color coded curves. The right 
panel shows the residuals between the observed and fit- 
ted color excess. There is a slight tail to the distribution 
owing to the lower quality of the fit to the high stellar 
mass, high metallicity and high color excess data. The 
rms dispersion between the observations and the fit is 
-0.11 dex. 

3 http://purl.com/net/mpfit 



We measure the SFRs for the DEEP2 sample of galax- 
ies using the H/3 luminosity. In order to obtain a flux 
calibration, we integrate the DEEP2 spectra over the I- 
band filter response and compare to the observed /-band 
magnitude. This gives us an estimate for the conversion 
between counts on the detector and flux given in ergs 
s . This procedure accounts for both throughput and 
slit losses assuming that the line-to-continuum emission 
ratio is constant inside and outside of the slit. This as- 
sumption is reasonable as slit losses are generally small 
because the vast majority of DEEP2 galaxies have effec- 
tive radius smaller than the 1" slit width (Weiner et al. 
2007). We linearly fit the flux calibrated continuum in 
a 80A window about the H/3 emission line and scale the 
equivalent width by the continuum to obtain a H/3 flux. 

We use the extinction law of Cardelli et al. (1989) and 
the fitted E(B— V) values given by Equation 6 to cor- 
rect the observed H/3 flux in the DEEP2 sample for dust 
attenuation. We obtain the luminosity from the flux us- 
ing the relation L(R(3) = F^p x ^D 2 L where F-^p is the 
dereddened H/3 flux and Dl is the luminosity distance 
determined from the redshift. We determine the SFR 
from the Kennicutt (1998a) relation between SFR and 
Ha luminosity given by 

SFR = 7.9 x 10" 42 i(Ha)[ergs s -1 ]. (7) 

We take L(Ha) = 2.86 x L(H/3) and scale the SFR from 
Salpeter to Chabrier IMF by dividing by 1.7. 

E06 infer their SFRs from the Ha luminosity corrected 
for dust extinction using the E(B— V) determined from 
SED fitting. A factor of 2 aperture correction was ap- 
plied and the SFRs were determined using the Kenni- 
cutt (1998a) relation converted to a Chabrier IMF. We 
redetermine the SFRs for their 87 galaxies by applying 
an extinction correction determined from Equation 6 to 
their observed Ha luminosities. Following E06, we apply 
a factor of 2 aperture correction to the observed lumi- 
nosity. We determine the SFR from Equation 7 and bin 
data according to E06. 

Figure 2c shows the distribution of the SFRs for all 
three samples. 

4. GALACTIC STELLAR MASS GROWTH 

4.1. The Observed Relation between Stellar Mass and 

SFR 

In Figure 4 we plot the MS relation for the SDSS (blue 
squares), DEEP2 (black circles) and E06 (red triangles) 
star-forming galaxies at redshifts of z = 0.07, 0.78 and 
2.26, respectively. For the SDSS, DEEP2 and E06 sam- 
ples the median SFR is determined in bins of stellar 
masses. For the SDSS, DEEP2 and E06 samples there 
are ~ 700, ~ 90 and — 15 galaxies in each data bin. 
The flattening in the relation observed at lower stellar 
masses is due to incompleteness (see below) such that 
the lower SFR galaxies fall out of our selected samples 
at lower stellar masses. Down to log(M») ~ 9.2 for the 
SDSS and DEEP2 sample and log(M*) ~ 9.5 for the E06 
sample the relation between stellar mass and SFR is well 
described by a power law. 

For the SDSS and DEEP2 samples the error bars plot- 
ted on the data points indicate the interval containing 
68% of the data within each mass bin. For the E06 sam- 
ple the error bars indicate the standard deviation of the 
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Fig. 4. — The logarithm of the SFR (Mq yr _1 ) plotted in bins 
of stellar mass for the SDSS (blue squares), DEEP2 (black circles) 
and E06 (red triangles) samples. For the SDSS and DEEP2 sample 
the median SFR is plotted and the error bars indicate the interval 
containing 68% of the data. For the E06 sample the median is 
plotted and the error bars are the standard deviation of the distri- 
bution. The dashed blue, solid black and dotted red lines are fits 
to the to the observed MS relation for the SDSS, DEEP2 and E06 
samples, respectively. The fit parameters are given in Table 2. 



TABLE 2 
Mass- SFR Relation 



Sample 


redshift 


5 


7 


SDSS 


0.07 


0.317 ± 0.003 


0.71 ± 0.01 


DEEP2 


0.78 


0.787 ± 0.009 


0.67 ± 0.02 


E06 


2.26 


1.608 ± 0.028 


0.46 ± 0.07 



Note. — The fitted MS relation for our three samples 
as parameterized by Equation 8. 



distribution within each mass bin. The median scatter 
in the SFRs is 0.25, 0.27 and 0.24 dex for the SDSS, 
DEEP 2 and E06 samples, respectively. The scatter in 
the three samples is comparable and is roughly constant 
with stellar mass. The observed scatter in our samples is 
comparable to the < 0.3 dex found in studies examining 
the M*-SFR relation out to z ~ 2 (Salim et al. 2007; 
Elbaz et al. 2007; Noeske et al. 2007b; Daddi et al. 2007; 
Pannella et al. 2009). 

We fit a linear relation between log('I') and 
log(M,i»/MQ) in the mass range where the samples are 
monotonically declining using the Unfit. pro routine in 
IDL. The form of the fit is given by 

log(SFR) = S + 7 ■ [log(M*/M Q ) - 10], (8) 

where 5 is the logarithm of the SFR at 10 10 Mq and 7 is 
the power law index. The dashed blue, solid black and 
dotted red lines in Figure 4 are the fits to the SDSS. 
DEEP2 and E06 samples, respectively. We give the fit 
parameters and the errors in Table 2. The errors for the 
SDSS and DEEP2 fit parameters are bootstrapped. The 
errors for the E06 sample are determined from propagat- 
ing the standard error of each data bin, a/y/n, where 
a is the standard deviation of the distribution and n is 
the number of objects in each bin. The power law index 
of the relation is consistent (within the small errors) for 
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Fig. 5.— The SFR (Mq yr" 1 ) plotted against stellar mass for our 
metallicity selected SDSS sample (black points) along with larger 
comparison sample from SDSS. The black points and error bars are 
the observed MS relation and its scatter for our metallicity selected 
SDSS sample. The dashed black line is our fitted relation given by 
Equation 8 with parameters given in Table 2. The solid black curve 
is the MS relation for our larger comparison sample and the gray 
lines are the 68 and 95% contours. 

the SDSS and DEEP2 samples. We have applied a fac- 
tor of 2 aperture corrections to all galaxies in the E06 
sample. This correction is likely overcorrecting the lower 
mass galaxies which could explain the shallower slope 
observed in the MS relation. The zero point of the MS 
relation shows evolution to higher values with increasing 
redshift. 

4.2. Sample Incompleteness 

In all three of our samples the MS relation flattens out 
at lower stellar masses (see Figure 4). The galaxies in 
all three of the samples are selected such that chemi- 
cal properties can be determined from the spectra. De- 
termining chemical properties from the emission lines of 
star- forming galaxies generally requires high S/N obser- 
vations and therefore samples with determined metallic- 
ities are more susceptible to incompleteness. We demon- 
strate that the observed flattening is due to incomplete- 
ness by comparing our sample of ^20,000 galaxies for 
which we have well determined metallicities with a larger, 
more complete sample of SDSS star-forming galaxies. 

We select the larger comparison sample by first re- 
moving AGN using the [OIIIJ/H/3 vs. [Nil] /Ha diagram 
(Kauffmann et al. 2003; Kewley et al. 2006) and imposing 
the similar redshift restrictions (z < 0.1) as our metallic- 
ity selected SDSS sample. This gives us a parent sample 
^200, 000 star-forming galaxies. The aperture corrected 
SFRs given in the SDSS DR7 rely on a Bayesian tech- 
nique principally fitting Ha and H/3 lines to determine ex- 
tinction corrected SFRs (Brinchmann et al. 2004; Salim 
et al. 2007). We select galaxies from the parent sample 
by requiring a S/N of 5 in the Ha and H/3 emission lines. 
This gives us a sample of ^140,000 star-forming galaxies 
with z < 0.1 and detections of Ha and Hf3. We note that 
the distribution of SFRs in our comparison sample are 
largely independent of choice of S/N cuts and therefore 
the observed incompleteness in our metallicity selected 
SDSS sample likely results from our limit on the equiva- 
lent width of K/3. 
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Figure 5 shows the SFRs for both our metallicity se- 
lected SDSS sample of galaxies and the larger compari- 
son sample. The black data points with error bars and 
dashed line are the same as in Figure 4. The solid black 
line is the median SFRs of our comparison sample sorted 
into 75 bins of stellar mass and the gray curves are the 
68 and 95% contours of the data. The MS relation in the 
larger comparison sample has a slope of 0.65 which is 
sightly more shallow than our metallicity selected SDSS 
sample. This slope is consistent with the slope deter- 
mined by Salim et al. (2007) for MS relation of SDSS 
galaxies. The zero point of the larger sample is slightly 
lower as compared to our metallicity selected SDSS sam- 
ple. Over the range where we fit the MS relation, the 
difference in the SFRs for the two samples is < 0.1 dex 
and both relations show similar scatter. 

From this comparison it is clear that the flattening 
observed in the MS relation in Figure 4 is due to in- 
completeness. The slightly higher zero point of our 
metallicity selected SDSS sample suggests that there is a 
slight bias against low star-forming galaxies over all mass 
ranges. This is consistent with our expectation that the 
higher S/N and equivalent width of H/3 required in the 
spectra for chemical analysis selects against galaxies with 
weak lines and low levels of star formation. Despite this 
selection effect, Figure 5 demonstrates that incomplete- 
ness does not significantly bias our determination of the 
MS relation. 

Noeske et al. (2007b) examine the MS relation for 
DEEP2 galaxies in the redshift range of 0.2 < z < 1. 
They determine SFRs from the 24 fim flux and add in the 
SFRs determined from the DEEP2 emission lines with 
no extinction correction to account for unobscured star 
formation. They determine the fit to the MS relation in 
the range where their sample is > 95% complete. The MS 
relation determined by Noeske et al. (2007b, but see Dut- 
ton et al. (2010)) at z = 0.78 is consistent with our de- 
termination for our DEEP2 sample at the same redshift, 
despite the fact that our sample is selected differently 
and we have determined SFRs from extinction corrected 
emission lines. Our correction for dust extinction relies 
on the correlation between extinction, stellar mass and 
metallicity (Section 4.3). This necessarily introduces a 
correlation between stellar mass and SFR. The fact that 
the MS determined from our DEEP2 sample is consis- 
tent within the errors to the relation fit by Noeske et al. 
(2007b) implies that our dust correction is not introduc- 
ing a spurious correlation between stellar mass and SFR. 
Furthermore, we conclude that the MS relation we deter- 
mine at z = 0.78 from our DEEP 2 sample is not biased 
due to incompleteness. 

The completeness of the E06 sample at z = 2.26 is 
much more uncertain. This sample is UV-selected which 
may bias against star-forming galaxies with low levels of 
star formation if they are present at this redshift. How- 
ever, we note that the observed relation determined from 
the E06 sample is consistent within the errors to the zero 
point and scatter found by Daddi et al. (2007) using a 
highly complete sample at z ~ 2 and Pannella et al. 
(2009) who use radio stacking to determine SFRs. This 
suggests that the observed MS relation in the E06 sample 
is not significantly biased. 
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Fig. 6.— The SFR (Mq yr" 1 ) for a 10 10 Af Q galaxy (0 O , top 
panel) and the power law index of the MS relation (7, bottom 
panel) as a function of redshift for our data (black points) alongside 
several values from the literature. The data come from Noeske et al. 
(2007b, red triangles), Elbaz et al. (2007, blue squares), Salim et al. 
(2007, yellow star), Daddi ct al. (2007, cyan triangle) and Pannella 
et al. (2009, green diamond). The dashed line in the top panel is 
the fit for ipo as a function of redshift for our three samples given by 
Equation 11. The dashed line in the bottom panel is our adopted 
value of 0.7 for the power law index of the MS relation, 7. 



4.3. Parameterization of SFR 

With a well defined MS relation at several redshifts, 
we can parameterize the SFR as a function of mass and 
redshift. We seek a parameterization where the variables 
of mass and redshift are separated such that the SFR is 
given by 

*(M»,z) = ^ («)-/(M«). (9) 

Here \&(M*, z) is the SFR as a function of stellar mass 
and redshift, ip (z) is the zero point which evolves with 
redshift and /(M*) is the relation between stellar mass 
and SFR. This definition explicitly assumes that the 
slope of the relation between stellar mass and SFR has 
no time dependency. 

As discussed in Section 4.1, the SFR as a function of 
stellar mass is well fit by a power law in stellar mass 
with an index near unity. We therefore parameterize the 
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stellar mass dependency as 



f(AQ = (M*/10 1U M Q 



(10) 



where the only free parameter, 7, is the power law index 
of the relation. With this definition for /(M*), the zero 
point, ip {z), is the SFR of a 10 10 M Q galaxy. The zero 
point of the relation as a function of redshift, as we show 
below, is well fit by an exponential and we parameterize 
it as 

ipo(z) = a ■ exp(/3z). (11) 

Here, a is the zero point at z = and /3 determines the 
rate of exponential increase with redshift. 

Figure 6 shows a compilation taken from the literature 
for the values of 7 (bottom panel) and tp Q (top panel) as 
defined in Equations 10 and 11, respectively. For com- 
parison, data from this study are plotted as black points 
on these plots. We have converted the stellar masses and 
SFRs for the data shown in Figure 6 such that they are 
consistent with the (Chabrier 2003) IMF. The data selec- 
tion and method of SFR determination vary. In the local 
universe, Salim et al. (2007) derive SFRs by fitting the 
UV/optical broad-band SED. Elbaz et al. (2007) adopt 
the Brinchmann et al. (2004) aperture corrected SFRs 
given in the DR4. At 0.2 < z < 1, the measurements 
of Noeske et al. (2007b) 4 have SFRs determined from 
a combination of uncorrected emission line fluxes and 
24 /xm fluxes. At z ~ 1, Elbaz et al. (2007) determine 
SFRs from the infrared luminosity inferred from the 24 
/itm flux. At z ~ 2, Daddi et al. (2007) determine SFRs 
from dust corrected UV luminosities and Pannella et al. 
(2009) determine SFRs from stacking analysis of the ra- 
dio continuum. 

Despite the differences in methodology, our determi- 
nation of the zero point, tp (top panel Figure 6), is re- 
markably consistent with values found in the literature. 
In the local universe, we find a marginally higher zero 
point, which is due to incompleteness (see Section 4.2). 
At z = 0.78 we determine the same relation (within the 
errors) as Noeske et al. (2007b). At z — 2.26 the rela- 
tion determined from the E06 sample is consistent with 
the higher redshift data. Out to z ~ 1, the slope of the 
MS relation, 7 (bottom panel of Figure 6), are consistent 
in most of the studies. At z > 1 the slopes show some 
variation, though the variation is only about ±35%. We 
note that in this study, our results are not sensitive to 
the adopted value of 7 because the redshift evolution of 
the SFR is significantly greater than the magnitude of 
the variation introduced by a ±35% change in 7. 

In order to derive ^(M^,z), we need to determine 
values for the three free parameters, a,/3 and 7. For 
the SDSS and DEEP2 samples the slopes are consistent 
(within the errors) and we adopt a constant value of 0.7 
for 7 in Equation 10 (see also Whitaker et al. 2012). This 
is shown by the dashed line in the bottom panel of Figure 
6. To determine a and B we fit a line of the form 



log[^o(*)] =log(a) + 0.434 pz, 



(12) 



to our data. The fit is shown by the dashed line in the 
top panel of Figure 6. We note that we have fit only to 
our data, but the fit is more or less consistent with all 

4 We have taken the parameters of the Noeske et al. (2007b) MS 
relation from the compilation of Dutton et al. (2010). 
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Fig. 7. — The fitted relation for the SFR (Mq yr' 1 ) as function 
of redshift, >&(M*, z), given by Equation 13 plotted over our three 
samples. 

the studies compiled from the literature. At z > 2.26 
we assume that ip a is flat noting that our results are not 
strongly dependent on this extrapolation because in our 
model most of the stellar mass growth in star-forming 
galaxies in the local universe occurs at z < 2.26 (see also 
Leitner 2012). 

With both ijj {z) and f(M*) determined, we can pa- 
rameterize the SFRs of star-forming galaxies as a func- 
tion of stellar mass and redshift by entering the fitted 
results back into Equation 9. The SFR as a function of 
stellar mass and redshift is given by 



*(Af»,z) = 2.00-exp(1.33z) 



10 



0.7 



,10 



[Meyr- 1 ]. (13) 



In Figure 7 we show the fitted relation plotted over the 
observations. Though the fitted slope of the MS rela- 
tion at z = 2.26 was shallower (see Table 2), the data 
is still reasonably well described by the relation given in 
Equation 13. 

5. INVENTORY OF OXYGEN 

In this study we compare an estimate of the total 
amount of oxygen produced with the amount found in 
the ISM and stars of local star-forming galaxies. We fo- 
cus on oxygen because it is a primary element produced 
in massive stars that end their lives as Type II super- 
novae. Oxygen is the most abundant element produced 
and therefore the easiest to measure from strong nebular 
emission lines originating in the ISM of galaxies. Oxygen 
production depends only weakly on the initial metallic- 
ity (Thomas et al. 1998; Kobayashi et al. 2006) and is 
considered a robust tracer of total metal production. We 
assume that oxygen production is independent of metal- 
licity. The rate of oxygen production is then simply given 
by 

dMZ 



dt 



(14) 



Here M£ is the mass of oxygen produced in solar mass 
units, P Q is a constant representing the mass of newly 
synthesized oxygen per solar mass of gas converted to 
stars and VP is the SFR. We can estimate the total oxygen 
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production by integrating Equation 14. This is given by 



M°(t > ti) =P J ^(t')dt'. 



(15) 



The amount of oxygen produced between some initial 
time, ti, and some later time, t, is proportional to the 
total amount of star formation during that epoch. This 
equation implies that M£ oc M*, where M£ is the to- 
tal amount of oxygen produced in a galaxy and M„ is 
the stellar mass of the galaxy. ^ is explicitly parame- 
terized as a function of redshift (see Equation 13) which 
we convert into a function time by assuming the con- 
version between redshift and time using the IDL routine 
galage.pro. 

The second component of our oxygen census is the oxy- 
gen found in the ISM. We define the global gas-phase 
oxygen mass abundance such that Z g = M°/M g . Here, 
M g and M° are the mass of gas and the mass of oxy- 
gen in the gas-phase, respectively, given in solar mass 
units. The total oxygen mass in the gas-phase is given 
by rewriting this equation such that M° = Z g M g . Both 
Z g and M g are measured quantities. 

The final component of our census of oxygen in star- 
forming galaxies is the oxygen locked up in stars. Assum- 
ing that the gas in galaxies is well mixed and mass return 
from stars is instantaneous, the rate at which oxygen is 
locked up in stars is given by 



dMl 
dt 



(1 



R)Z g V. 



(16) 



Here M° is the mass of oxygen going into stars in solar 
mass units and R is a constant representing the fraction 
of gas that is returned to the ISM through various stellar 
mass loss processes. We can integrate Equation 16 to get 
the total amount of oxygen locked up in low mass stars. 
This is given by 

M,°(i > U) = (1 - R) j Z g (t') *{t') dt'. (17) 

M° represents the amount of oxygen in the ISM that gets 
forever locked up in low mass stars in the time between 
ti and t. Similar to Z g is a explicitly a function of 
redshift which we convert into a function of time. In 
Sections 6 and 7 we develop an empirically constrained, 
self-consistent approach for determining ^(t) and Z g (t), 
respectively. 

6. TRACING GALAXIES THROUGH COSMIC TIME 

One interpretation of the MS relation and its evolution 
is that most star formation is driven largely by secular 
processes such as gas accretion and mergers do not play 
a significant role (e.g. Noeske et al. 2007b; Dutton et al. 
2010; Elbaz et al. 2011). Under the assumption that 
stellar mass build up is a purely secular process, we can 
trace stellar mass evolution of galaxies by assuming that 



dM* 
dt 



(18) 



Here is the time rate of change of stellar mass, R is 
the return fraction which accounts for gas that is formed 
into stars but then returned to the ISM via various mass 
loss processes (e.g. supernovae and stellar winds) and \I/ 



is the SFR. For simplicity we assume that R is a constant 
(see Section 9 and the Appendix for more discussion). 

If we further assume the continuity condition that star- 
forming galaxies evolve along the MS relation and SFRs 
are a (determined) function of stellar mass and redshift 
(see Equation 13), it is straightforward to integrate Equa- 
tion 18 in order to determine the stellar mass as a func- 
tion of redshift. This is given by 



M,(t) = M.,i + (l 



- R) i 



V(M*,t')dt'. (19) 



Here M*(t) is the stellar mass at some redshift, t < ti. 
M* t i is the stellar mass at some initial time, ti, and the 
integral of ^(M^,t) represents the gain in stellar mass, 
modulo the return fraction R, from some initial time ti 
to some later time t. If M*^ is set to an arbitrarily low 
value (10 6 M Q ) and the upper limit of the integral is the 
current epoch, then ti is simply the formation time. This 
is similar to the staged model of Noeske et al. (2007a) 
such that lower stellar mass star-forming galaxies in the 
local universe begin forming their stars at later times. 

Equation 19 is a simple, yet powerful tool for identi- 
fying star-forming galaxies at lower redshifts with their 
progenitors at higher redshifts under the assumption of 
secular evolution (i.e. no contribution from merging) and 
that the SFRs are a known function of stellar mass and 
redshift (see Equation 13); the latter being observation- 
ally determined. Similar approaches have been devel- 
oped in the literature and we direct the reader to Lcit- 
ner (2012) for a detailed discussion of the implications 
of such a model for stellar mass growth. Because of the 
stellar mass dependence of <J/ and the fact that SFR is 
measured in units of solar masses per year and not solar 
masses per redshift, Equation 18 does not have a straight- 
forward analytical solution and thus requires numerical 
integration. The integration of Equation 18 gives both 
the stellar mass and star formation history of the galaxy. 

7. GALACTIC CHEMICAL EVOLUTION 

The three data samples used in this study have been 
selected because they all have well determined MZ rela- 
tions. The MZ relation for the SDSS and DEEP2 sam- 
ples are determined by Zahid et al. (2011) and for the 
E06 sample by Erb et al. (2006a). We convert the metal- 
licity, traditionally quoted as a number density given by 
12 + log(0/H), to gas-phase mass abundance using the 
simple linear conversion given by Equation 6. 

In Figure 8 we plot the MZ relation for our three data 
samples. The MZ relation for the SDSS and DEEP2 sam- 
ple are determined using the calibration of Kobulnicky & 
Kewley (2004, Figure 8a) and the for the E06 sample us- 
ing the calibration of Pettini & Pagel (2004, Figure 8b). 
Kewley & Ellison (2008) have shown that while various 
calibrations give relatively accurate measurements of the 
mctallicity, there are systematic differences. We apply 
the conversions constants they derive in order to convert 
between the two calibrations shown in Figure 8. The 
shape of the MZ relation determined from the two cal- 
ibrations is very similar. The main difference being the 
0.35 dex higher abundances of the Kobulnicky & Kewley 
(2004) calibration with respect to the calibration of Pet- 
tini & Pagel (2004). In this study we adopt the average 
of the two calibrations for our determination of metallici- 
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Fig. 8. — The MZ relation determined using the calibration of a) Kobulnicky k, Kewley (2004) and b) Pettini & Pagel (2004). The 
metallicity is given as the mass abundance of oxygen relative to hydrogen. The blue, black and red curves are the MZ relations determined 
for the local sample from SDSS, the intermediate redshift sample from DEEP2 (Zahid et al. 2011) and the high redshift sample from Erb 
et al. (2006a), respectively. The metallicities for the SDSS and DEEP2 data are originally determined using the Kobulnicky & Kewley 
(2004) diagnostic and the metallicity of the Erb et al. (2006a) sample is originally determined using the Pettini & Pagel (2004) diagnostic. 
When necessary, the metallicities have been converted using the coefficients given by Kewley & Ellison (2008). 



ties. In Section 8.3 we provide a more detailed discussion 
of the calibrations uncertainties. 

We use the observed MZ relations to determine the gas- 
phase mass abundance as a function of stellar mass and 
redshift. The three samples in this study suggest that 
the MZ relation evolves linearly with time. Moustakas 
et al. (2011) investigate the MZ relation evolution out to 
z = 0.75 and suggest that it evolves linearly with red- 
shift. In this study we assume linear evolution with time 
but note that assuming linear evolution with redshift has 
only minor quantitive effects on our results and has no ef- 
fect on our interpretation. We further assume a constant 
offset between the three relations at log(Af*/M Q ) < 9.2. 
The MZ relation in the local universe is shown to extend 
to low stellar masses (Lee et al. 2006; Zahid et al. 2012) 
and we linearly extrapolate the three relations to lower 
stellar masses when necessary. 

In summary, we determine the gas-phase mass abun- 
dance of oxygen at arbitrary redshifts and stellar masses 
by linearly interpolating between the three observed rela- 
tions shown in Figure 8. We linearly interpolate both in 
time and stellar mass in order to determine Z g (M*, z). 

8. A CENSUS OF OXYGEN 

There are three components to our oxygen census: to- 
tal oxygen produced, amount of oxygen locked up in stars 
and the amount of oxygen in the gas-phase. To put this 
in context, we begin our census by noting the basic bal- 
ance equation for oxygen. In a closed-box system with 
no inflows or outflows of gas, the oxygen balance is given 
by 



RZg^ 



= P *-(l-R)Z g V 



(20) 



As before, M° is the mass of oxygen in the gas-phase, 
Z g is the gas-phase mass abundance of oxygen and is de- 
fined as Z g = M°/Mg, where M g is the gas mass. \& 
is the SFR. R and P are the gas return fraction and 
mass of newly synthesized oxygen per unit stellar mass, 



respectively. The first term on the right-hand side of 
Equation 20 represents the total mass of oxygen pro- 
duced, the second term represents the amount that goes 
into stars and the third term represents the amount that 
is subsequently returned through mass loss. The stel- 
lar yield, y, is related to the nucleosynthetic yield and 
gas fraction by y = P a /(l — R)- We introduce this ba- 
sic balance equation to place the census of oxygen in 
star-forming galaxies that follows within the context of 
simple chemical evolution models. For a more detailed 
discussion of these chemical evolution models we refer 
the reader to Edmunds (1990). 

8.1. Total Oxygen Production 

The total amount of oxygen produced by the stars 
within a galaxy is simply given by 



Po 



l-R 



(21) 



The quantity M*/(l — R) represents the total amount 
of gas converted into stars. Several simplifying assump- 
tions go into Equation 21: 1) The IMF is constant, 2) 
gas is instantaneously returned and 3) the oxygen yield 
is independent of metallicity. In Section 8 we discuss sys- 
tematics and uncertainties associated with each of these 
assumptions. 

8.2. Oxygen in the Gas-Phase 

The total mass of oxygen in the gas-phase in local star- 
forming galaxies is given by M° = Z g M g , where M° 
is the mass of oxygen in the gas-phase. Unfortunately, 
direct gas mass measurements are available for only a 
few galaxies in our local sample. We therefore estimate 
the gas mass of local galaxies using the observed relation 
between gas fraction and stellar mass for star-forming 
galaxies. Peeples & Shankar (2011) compile gas mass 
estimates (HI + H2, with correction for helium) from 
several sources in the literature and give the binned cold 
gas fraction as a function of stellar mass (see Table 2 
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Fig. 9. — The gas fraction as a function of stellar mass for star- 
forming galaxies taken from the compilation of Peeples & Shankar 
(2011). The dashed line is a fit to the relation and is given by 
Equation 22. 

in their paper). The values are compiled from the data 
sets of McGaugh (2005), Leroy et al. (2008) and Garcia- 
Appadoo et al. (2009). 

In Figure 9 we plot the log of the gas fraction, f g , as a 
function of stellar mass. The gas fraction is defined such 
that f g = Mg/M* 5 . The stellar mass estimates of Peeples 
& Shankar (2011) arc taken from MPA/JHU catalog and 
we make a 0.2 dex correction for consistency (Zahid et al. 
2011). The error bars are the la uncertainties on the 
mean of log(/ s ). The dashed line is a fit to the relation 
given by 

!og(/ fl ) = (-O.28±O.O3)-(O.55±O.O3)-[log(M,/M )-lO]. 

(22) 

The errors in the parameters come from propagating the 
uncertainties in the mean of log(/ 9 ). When calculating 
gas mass, we add 0.1 dex to the fitted relation at all stel- 
lar masses to account for ionized gas (see Table 1.1 in 
Tielens 2005). Using this relation, the gas mass is given 
by M g = f g M*. We note that the gas fractions given 
by Equation 22 are greater than values inferred from in- 
version of the Schmidt-Kennicutt relation (see Figure 2 
in Peeples & Shankar 2011). Using gas fractions de- 
termined from inverting the Schmidt-Kennicutt relation 
would therefore lower the estimate of the mass of oxygen 
in the gas-phase and increase the oxygen deficit presented 
in Section 8.4. 

We can straightforwardly estimate the mass of oxy- 
gen in the gas-phase of local star-forming galaxies from 
observed quantities by combining the relation for f g 
given in Equation 22 with the MZ relation for local star- 
forming galaxies shown in Figure 8. This is given by 

mi = y.,j,\\.. 

8.3. Oxygen Locked Up in Stars 

The most difficult component of our oxygen census to 
estimate is the oxygen locked up in stars. The total 
amount of oxygen produced by stars and the amount of 
oxygen found in the gas-phase of galaxies can be inferred 

5 We have adopted the convention of Peoples & Shankar (2011) 
in defining f g . This differs from the traditional definition of the 
gas fraction given by f g = M g /(M„ + M g ). 
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Fig. 10. — The a) stellar mass, b) SFR (M Q yr _1 ), c) gas-phase 
oxygen abundance and d) the mass of oxygen locked up in stars 
determined from our models as a function of look back time and 
redshift. The solid, dashed and dotted curves are galaxies that 
have a stellar mass of 10 11 , 10 10 and 1O 9 M in the local universe 
(z = 0.07). 

from the present-day physical quantities of the stellar 
mass and gas-phase abundance. However, the oxygen 
locked up in stars is an accumulated property that is de- 
pendent on the star formation and chemical history. The 
amount of oxygen locked up in stars between some initial 
redshift ti and some later redshift t is analytically given 

by 

M°(t > U) = (1 - R) J Zg(t') f(t') dt'. (23) 

We have parameterized both the SFR (see Section 4) 
and gas-phase oxygen abundance (see Section 7) as a 
function of stellar mass and time and are therefore able to 
empirically constrain the total amount of oxygen locked 
up in stars by simply carrying out a numerical integration 
of Equation 23. 

We simultaneously determine the star formation and 
stellar mass history of a galaxy using Equation 19. We 
set M t> i — 1O 6 M0 and carry the integral from ti out 
to t = 12.5 Gyr 6 . Higher values of ti correspond to 
younger, lower stellar mass star-forming galaxies in the 
local universe. We show the stellar mass history of three 
galaxies (ti = 2.4, 4.0 and 6.4 Gyr) in Figure 10a and 
star formation history in Figure 10b assuming R = 0.35. 
The stellar mass history as a function of redshift allows 
us to determine the metallicity as a function of redshift 
for each galaxy using the procedure described in Section 
7. We plot the chemical history in Figure 10c. Entering 
the star formation history, <3>(t) (Figure 10b), and the 
chemical history, Z g (t) (Figure 10c), into Equation 23, 
we can determine the mass of oxygen locked up in stars 
as a function of redshift. This is shown in Figure lOd. 

Gallazzi et al. (2005) and Panter et al. (2008) examine 

6 We set the upper limit of the integral to a time corresponding 
to a redshift of z = 0.07 rather than 2 = because this is the 
median redshift of our local sample from SDSS. 
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Fig. 11. — The stellar mass— stellar metallicity relation for lo- 
cal star-forming galaxies. The solid black curve is the fitted rela- 
tion determined by Panter et al. (2008) and the dotted curves are 
the 16 and 84% contours of the distribution. The red dotted and 
green dot-dashed curves plot the stellar metal content of local star- 
forming galaxies calculated from Equation 24 using the Kobulnicky 
& Kewley (2004) and Pettini & Pagel (2004) abundance calibra- 
tions, respectively. For reference (c.f. Panter et al. 2008, Figure 6), 
the solid yellow curve plots the MZ relation from Tremonti et al. 
(2004). 

the stellar mass— stellar metallicity relation determined 
from stellar population modeling of SDSS spectra. Both 
groups derive a consistent relation which implies that the 
imprint of the MZ relation from earlier epochs is observed 
in stellar populations of galaxies today. In Figure 11 wc 
compare the metallicity of stars from our model analysis 
to that observed by Panter et al. (2008) using ~ 300, 000 
galaxies in SDSS. We estimate the mass-weighted stel- 
lar metallicity from our model of the star formation and 
chemical histories of galaxies (see Figure 10b and c, re- 
spectively). In particular we compute 

\l y(t')z(t')dt' r* y(t')z(t')dt' 

Z*{M*) = H ' t = U ' , (24) 

where is the mass-weighted stellar metallicity as a 
function of current stellar mass, M*. We compare di- 
rectly to the mass-weighted stellar metallicity relation 
observed and parameterized by Panter et al. (2008, see 
Equation 1 and Figure 6). In Figure 11 the Panter et al. 
(2008) fitted relation is plotted by the solid black curve 
and the 16 and 84% contours of the distribution are 
shown by the dotted black curve. 

As can be seen from Equation 24, our estimate of the 
metal content of stars is dependent on the choice of gas- 
phase abundance diagnostic. In Figure 11 we plot the 
mass- weighted stellar metallicity for both the Kobulnicky 
& Kewley (2004, red dashed curve) and Pettini & Pagel 
(2004, green dot-dashed curve) calibration. For reference 
the MZ relation from Tremonti et al. (2004) is plotted 
by the yellow curve and for consistency we adopt the so- 
lar metallicity value of 12+log(O/H) = 8.87 used by 
Panter et al. (2008). Estimates from the two different 
abundance diagnostics converted into stellar metallici- 
ties using Equation 24 envelope the relation from Panter 
et al. (2008). In this study we have adopted a metallicity 
calibration that is an average of the Kobulnicky & Kew- 



ley (2004) and Pettini & Pagel (2004) calibration (see 
Section 7). Within the systematic uncertainties, our es- 
timate of stellar metallicities using Equation 24 and the 
observed stellar metallicities from Panter et al. (2008) 
are consistent. 

8.4. The Oxygen Deficit 

In Figure 12a-c we plot the mass of oxygen for our 
three census components. We plot the total mass of oxy- 
gen produced (blue curve), the mass of oxygen in the 
gas-phase (black curve) and the mass of oxygen locked 
up in stars (red curve) as a function of stellar mass for 
star-forming galaxies in the local universe. The solid, 
dotted and dashed curves are used to illustrate the ef- 
fect varying the parameters has on different components 
of the census. The three free parameters in our model 
are the return fraction, R, the oxygen yield, P Q , and a 
constant offset in metallicity (which we refer to as AZ g ) 
since the zero point varies depending on the calibration. 
Figure 12a-c demonstrates the dependencies of the vari- 
ous components of our census on the free parameters of 
the model. In Figure 12a-c we vary R, P a and AZ g , re- 
spectively, while keeping the other two parameters fixed. 
The three values adopted for the varying parameter and 
the adopted values of the constant parameters are shown 
in a-c. In Figure 12a we see that variations in R mainly 
affects the total amount of oxygen produced with a small 
effect on the amount of oxygen locked up in stars. Figure 
12b shows that varying P Q only affects the total amount 
of oxygen produced and Figure 12c demonstrates that 
varying AZ g affects both the stellar and gas-phase com- 
ponents of the census, but has no effect on the total 
amount of oxygen produced. 

For a closed-box model, the total mass of oxygen pro- 
duced equals the mass of oxygen found in the gas-phase 
and stars. However if metals are lost from the system, 
the equality will not hold. Thus we have 

AM" = M£ - M° - M°, (25) 

where, as before, Mj, is the total mass of oxygen pro- 
duced and M° and M° is the mass found in the gas- 
phase and stars, respectively. AM", which we refer to as 
the oxygen deficit, represents the total amount of oxygen 
produced that is unaccounted for by the gas-phase and 
stellar components. In Figure 12d we plot the oxygen 
deficit which is the same for the three cases shown in Fig- 
ure 12a-c. The variation in the parameters are degenerate 
such that the oxygen deficit depends only on the variation 
in the quantity AZ g (l — R)/P , effectively reducing our 
three free parameters to one free parameter. We note that 
the stellar yield, y, is given by y = P /(l — R)- There- 
fore, our model is only sensitive to the changes in the 
zero point of the metallicity relative to the stellar yield. 
Higher (lower) values of AZ g /y correspond to a smaller 
(larger) oxygen deficit. To first order, this is true for all 
models of chemical evolution relying on these quantities. 

If we assume that all oxygen in galaxies originates from 
the stars within the galaxy, the oxygen deficit cannot be 
negative because a negative value would imply that there 
is more oxygen in the gas and stars than was produced 
by the galaxy. A physical model of the oxygen deficit 
requires a combination of the free parameters (AZ g (l — 
R)/P ) that yield AM° > 0. The dotted line in Figure 
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Fig. 12. — The census model of oxygen determined by a) varying the return fraction, b) the oxygen yield and c) the zero point of the 
oxygen abundance calibration while keeping the other two parameters constant, respectively, a-c) The total oxygen produced (blue line), 
oxygen found in the gas (black) and stellar phase (red) are plotted as a function of stellar mass for local star-forming galaxies. In a-c) 
we display the three values adopted for the varying parameter along with the constant values adopted for the other two parameters. The 
solid, dotted and dashed lines indicate how the mass of oxygen varies for the three different census components when adopting the three 
different values for the varying parameter, d) The oxygen deficit which is defined as the total oxygen produced minus the oxygen found in 
the gas and stellar phase. The oxygen deficit is the same for the oxygen masses determined in a-c). 



12d gives a negative oxygen deficit at Af* < 10 10 M Q , 
implying this is an unphysical model. 

The functional form of the oxygen deficit with respect 
to stellar mass varies according to AZ g (l — R)/P . How- 
ever, a generic feature of the oxygen deficit is that it 
monotonically increases with stellar mass and is substan- 
tial at high stellar masses regardless of our particular 
choice of AZ g (l — R)/P . In the case where the oxygen 
deficit is substantial at all stellar masses (solid curve in 
Figure 12), we have AM° cx M" with a near unity value 
of for the exponent (a = 1.14). When the oxygen deficit 
is large at all stellar masses (solid curve in Figure 12), 
the oxygen deficit does not strongly depend on the shape 
of the MZ relation or the shape of the relation between 
gas fraction and stellar mass. At more moderate values 
of AZ g (l — R) jP Q examined in our models (dashed curve 
in Figure 12), the relation between the oxygen deficit and 
stellar mass is more complicated. In Section 9.1 we argue 
that a model where the oxygen deficit is large at all stel- 
lar masses and is nearly proportional the stellar mass is 
most consistent with independent observations of oxygen 
in the halos of star-forming galaxies. 

In Figure 13 we plot the fractional oxygen deficit. In 
the top panel we plot the logarithm of 8 , which we define 
as 



MS 



M° 
g 



S a is the oxygen deficit relative to the total oxygen con- 
tent of the galaxy. In the bottom panel we plot the log- 
arithm of p , which we define as 



Po = AM°/M° 



(27) 



S a is an important quantity reflecting key physical pro- 
cesses governing chemical evolution of galaxies but can- 
not be determined for individual galaxies as the fraction 
of oxygen locked up stars cannot be measured. p a does 
not account for oxygen locked up in stars, but can be de- 
termined from observable quantities and is therefore an 
important parameter for comparison with observations. 

9. SYSTEMATICS AND UNCERTAINTIES 

We have derived the oxygen deficit by comparing the 
total oxygen produced with the amount found in the gas- 
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Fig. 13. — The top panel shows the oxygen deficit relative to the 
mass of oxygen in the gas and stellar phase and the bottom panel 
shows the oxygen deficit relative to only the mass of oxygen in the 
gas-phase. The oxygen deficit is the same as in Figure 12d and the 
solid, dotted and dashed lines indicate the different oxygen deficit 
derived by varying the free parameters. 

phase and stars of galaxies. We have derived the oxygen 
deficit on the basis of observed relations that are deter- 
mined from large samples of galaxies. Star formation 
rates and metallicities of individual galaxies are subject 
to stochastic variations resulting from physical processes 
that are not fully understood. The results of our study 
are valid in a statistical sense and do not necessarily ap- 
ply to any individual galaxy. 
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Large stochastic variations in the star formation rates 
and metallicities of galaxies does pose a problem for our 
approach since we assume galaxies evolve along the ob- 
served MS and MZ relations. The tight relation between 
stellar mass and SFR and the small scatter suggest that 
the SFRs of star-forming galaxies in the local universe 
do not vary rapidly over their lifetime and that mergers 
play a minor role (e.g. Noeske et al. 2007b; Dutton et al. 
2010; Elbaz et al. 2011). Though the observed scatter in 
metallicities increases at the low stellar mass end of the 
local MZ relation (Zahid et al. 2012), the la scatter is 
< 0.2 dcx (Trcmonti et al. 2004). The main source of un- 
certainty then is likely not associated with our approach 
but rather with parameters of our model and other sim- 
plifying assumptions which we discuss below. 

9.1. The Return Fraction and Instantaneous Recycling 
Approximation 

The fraction of gas returned to the ISM is mass depen- 
dent and therefore sensitive to the particular choice of an 
IMF, with theoretical values ranging between ~ 20—50%. 
We have adopted this range in Figure 12a. The physical 
mechanism regulating mass loss in stars varies by stel- 
lar mass. In low mass stars, stellar winds during the 
red-giant branch and asymptotic giant-branch are the 
dominant mechanism. In intermediate mass stars winds 
during the asymptotic giant branch phase dominate and 
for high mass stars, both stellar winds and supernovae 
contribute to stellar mass loss (Jungwiert et al. 2001, and 
references therein). 

The return fraction is important in our models in 
that it sets the total amount of star formation required 
for a galaxy to have a particular current stellar mass. 
The total amount of gas converted into stars, M gs , in 
a galaxy with a current stellar mass of M* is given by 
M gs = M*/(l - R). A galaxy with R = 0.5 will, to first- 
order, form 40% more stars than a galaxy with R = 0.3 
to achieve the same stellar mass. Because the amount 
of oxygen produced is directly proportional to the total 
amount of star formation, there will be a commensurate 
increase in total oxygen production. 

We have also adopted the instantaneous recycling ap- 
proximation (IRA) introduced by Tinsley (1980). The 
IRA assumes that the gas is returned to the ISM imme- 
diately following star formation. It is shown that the IRA 
works well for high mass stars whose lifetimes are short 
compared to star formation time scales. Oxygen is a pri- 
mary element, mainly produced in Type II supernovae. 
In our models, the lifetimes of massive stars responsi- 
ble for oxygen enrichment and gas mixing timescales are 
significantly shorter than the timescales of star forma- 
tion and therefore the IRA is reasonable in regards to 
oxygen production. However, the gas return from Type 
II supernovae accounts for 10-20% of the gas returned 
to the ISM and a larger fraction is returned on longer 
timescales from intermediate and low mass stars (see 
Jungwiert et al. 2001). 

Examining various IMFs, Leitner & Kravtsov (2011) 
show that the majority of stellar mass loss occurs in the 
first 1-2 Gyrs (c.f. their Figure 1). Noeske et al. (2007b, 
among others) argue that the MS relation implies that for 
most star-forming galaxies the SFR gradually declines. 
The gradual decline in SFRs is a generic feature of mod- 
els assuming the continuity condition that galaxies build 



up stellar mass evolving along the observed MS relation 
(e.g. Peng et al. 2010; Leitner & Kravtsov 2011; Leitner 
2012, this work). The IRA can be reasonably applied in 
situations where the SFR is not changing rapidly with 
rcdshift because though the current generations of stars 
being formed will return majority of their gas in 1-2 Gyrs, 
the total gas return will be an integration over previous 
generations. If the SFR is not rapidly varying, the re- 
turn fraction will reach a relatively constant value within 
a few Gyrs. 

In the Appendix we apply a time dependent treatment 
for gas return in determining the oxygen deficit. We no 
longer assume the IRA so that oxygen is not instanta- 
neously recycled nor is the return fraction a constant. 
We show that our results are not sensitive to time de- 
pendencies in the return rate and adopting a constant 
and instantaneous return of gas is valid. 

9.2. The Oxygen Yield 

The production of oxygen in our models is defined by 
the oxygen yield, P Q . The oxygen yield represents the 
total amount of newly synthesized oxygen per unit stellar 
mass of star formation. Following Henry et al. (2000), we 
formally define the integrated oxygen yield as 

P a = mp {m)(j){m)dm. (28) 

J mi 

Here m; and m u are the lower and upper limits of the 
mass range of stars formed, respectively, p {m) is the 
mass fraction of oxygen synthesized by a star of stellar 
mass m and <j)(m) is the normalized IMF. Uncertainties 
in P Q arise from uncertainties in the stellar yields, the 
IMF and the mass range of integration. Here we discuss 
only the stellar oxygen yields. For a more complete re- 
view of uncertainties in stellar yields of various elements 
see Romano et al. (2010). 

The dependence of the oxygen yield on the IMF is clear 
from examination of Equation 28. The relative weight of 
stars in different mass ranges assuming different IMFs 
are given by Romano et al. (2005, Table 1). The rel- 
ative weight of 8 — 40M Q stars is a factor 1.6 higher 
for the Chabrier IMF as compared to the Salpeter IMF, 
highlighting the uncertainties in the oxygen yields asso- 
ciated with the assumption of the IMF. These uncertain- 
ties are further exacerbated by the possibility of a vary- 
ing IMF. Kroupa & Weidner (2003) argue that the inte- 
grated galactic initial mass function (IGIMF) , the stellar 
mass distribution function for all stars within a galaxy, 
is an integration of the stellar IMF over the embedded 
cluster mass function, ECMF. There exists an empirical 
relation between the ECMF and the observed SFR, lead- 
ing to a dependence of the IGIMF on the observed SFR. 
Peoples & Shankar (2011) show several models for the 
stellar yield, taking into consideration these variations. 

Almost all oxygen in the universe is produced in mas- 
sive stars (m > 8M ). Using detailed models of Type II 
supernovae explosions, Woosley & Weaver (1995) provide 
the nucleosynthetic yields for stars in the mass range of 
11 — 4OM and for a range of metallicities. Thielemann 
et al. (1996) present similar model calculations but in the 
mass range of 13 — 25M . Thomas et al. (1998) com- 
pare these two models finding good agreement between 
the oxygen yields except at high stellar masses where 
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the yields calculated by Woosley & Weaver (1995) sat- 
urate. Henry et al. (2000) calculate the integrated oxy- 
gen yield using the models of Maeder (1992), Woosley & 
Weaver (1995) and the Thielemann et al. (1996) mod- 
els updated to include contributions from 40 and 70M Q 
stars (Nomoto et al. 1997). Assuming a Salpeter IMF 
and a range of upper and lower masses Henry et al. (2000) 
find that the oxygen yields range between 0.004 — 0.016 
(see their Table 2). Using more sophisticated models, 
Kobayashi et al. (2006) calculate nucleosynthetic yields 
by assuming a range of metallicities, explosion energies 
(supcrnovae and hypernovae) and mctallicity dependent 
mass loss. They calculate the integrated oxygen yield 
using a Salpeter IMF and a lower and upper mass limit 
of 0.07 and 50M Q , respectively. They find that the oxy- 
gen yield is weakly dependent on metallicity and ranges 
between 0.006 — 0.008. In Figure 12b, we adopt a range 
of 0.0056 - 0.0091 for the oxygen yield. 

9.3. The Abundance Calibration 

The uncertainty in the absolute nebular abundance cal- 
ibration is a long-standing problem in observational as- 
tronomy. Kcwley & Ellison (2008) show that the metal- 
licity can vary by as much as 0.7 dex when using differ- 
ent abundance diagnostics for the same set of galaxies. 
This poses serious problems for understanding galactic 
chemical evolution. The uncertainties are largely due to 
method of calibration used in establishing the diagnos- 
tic. Empirical methods rely on calibrating strong line 
ratios using HII regions and galaxies with metallicities 
determined using the so-called direct method which uses 
temperature sensitive auroral lines, (e.g. Pettini & Pagel 
2004). Empirical methods have several limitations and 
theoretical calibrations of strong-line ratios using pho- 
toionization models have also been developed (e.g. Kew- 
ley & Dopita 2002; Kobulnicky & Kewley 2004; Tremonti 
et al. 2004) . Detailed discussions of strengths and weak- 
nesses of the various methods and calibrations can be 
found in Kewley & Ellison (2008) and Moustakas et al. 
(2010), we summarize the salient points here. 

Pcimbcrt (1967) points out that temperature fluc- 
tuations within nebular regions may lead to overesti- 
mates of the temperature. Temperature fluctuations are 
thought to be more problematic in metal-rich HII re- 
gions, where efficient line-cooling may lead to tempera- 
ture inhomogeneities and strong temperature gradients 
(Garnett 1992). These overestimates of nebular temper- 
atures lead to underestimates of the metallicity with the 
direct method (Stasifiska 2002, 2005; Bresolin et al. 2006) 
and temperature inhomogeneities among the HII regions 
within a galaxy also pose problems for metallicity de- 
terminations based on global spectra (Kobulnicky et al. 
1999). A second, and perhaps more pernicious prob- 
lem is the distribution of electrons within an HII region 
may not follow a Boltzmann distribution but rather may 
be K-distributed (a well known particle distribution in 
plasma physics). A K-distribution of electrons may ac- 
count for the discrepant temperatures and metallicities 
inferred from auroral lines (Nicholls et al. 2012). 

Empirical calibrations require observations of ex- 
tremely faint auroral lines. In particular, the [OIIIJA4363 
line is on the order of 100 times weaker than the strong 
optical oxygen lines (e.g. [OIIJA3727, [OIII] A5007). Fur- 
thermore, the line strength diminishes with increasing 



metallicity, due to efficient line cooling, and is only ob- 
served in low metallicity HII regions (< O.5Z ). The 
high S/N spectra required to observe the auroral lines 
taken together with temperature inhomogeneities within 
individual HII regions suggests that samples used to cal- 
ibrate empirical methods may be biased and unreliable, 
particularly at higher metallicities. 

Theoretical methods rely solely on photoionization 
models to calibrate strong-line ratios. These methods 
are not susceptible to observational limitations imposed 
by empirical calibrations and metallicities are well con- 
strained and parameterization is well defined over the full 
range of observed line ratios. The absolute calibration 
of metallicities are model dependent and uncertainties 
are subject to the simplifying assumptions of the model. 
Indirect evidence suggests that these methods system- 
atically overestimate the metallicity (e.g. Bresolin et al. 
2009a; Kudritzki et al. 2012). 

The metallicities in our samples are determined using 
cither the theoretical diagnostic of Kobulnicky & Kewley 
(2004, SDSS and DEEP2) or the empirical diagnostic of 
Pettini & Pagel (2004, E06). These two methods differ 
by a constant offset of <~ 0.3 dex. For our mctallicity 
diagnostic we therefore adopt the average of the metal- 
licities determined using these two methods. 

Kewley & Ellison (2008) have shown that the relative 
estimates of the metallicity are robust for various diag- 
nostics. However, the absolute calibration varies signifi- 
cantly and the shape of the MZ relation can change de- 
pending on the calibration adopted (see Figure 2 of Kew- 
ley & Ellison 2008). In our models, the oxygen deficit is 
sensitive to the shape of the MZ relation only in the case 
when the oxygen deficit is small. As we show in the fol- 
lowing section, our model showing a large oxygen deficit 
(the solid line in Figure 12d) is most consistent with ob- 
servations of oxygen in the halos of star-forming galaxies. 
Furthermore, the shapes of the MZ relations determined 
from the calibration of Kobulnicky & Kewley (2004) and 
Pettini & Pagel (2004) are consistent within the errors 
(see Figure 8) despite the fact that the two diagnostics 
are calibrated independently using a theoretical and em- 
pirical calibration, respectively. 

9.4. Depletion onto Dust Grains 

Some uncertainty in the absolute metallicity scale is 
associated with depletion of metals onto dust. Several 
authors have inferred ^0.1 dex depletion of oxygen onto 
dust grains (Jenkins 2004; Cartledge et al. 2006; Pcim- 
bert & Peimbert 2010). However, in NGC 300 the agree- 
ment between the gradient inferred from A and B super- 
giants and that from emission lines suggests that there is 
very little, if any dust depletion in this galaxy (Bresolin 
et al. 2009a). It is likely that dust depletion ranges be- 
tween ~ 0.1 dex and perhaps is dependent on physical 
properties of galaxies such as metallicity or stellar mass. 

The oxygen deficit we estimate is sensitive to the gas- 
phase abundance only at small oxygen deficits. If we 
assume a constant relative level of depletion (e.g. ^0.1 
dex), the effects of dust depletion on the oxygen deficit 
is more pronounced in lower mass galaxies which have a 
smaller oxygen deficit. However, low mass galaxies may 
also have lower levels of depletion as NGC300 7 indicates. 

7 The stellar mass of NGC300 is ~ 1O 9 M (Kent 1987). 
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By comparing to independent observations of the mass 
of oxygen found in the halos of star-forming galaxies, we 
argue in Section 10.1 that a model with a large oxygen 
deficit is favored. In the case of a large oxygen deficit 
(solid curve in Figure 12), assuming that oxygen is de- 
pleted onto dust by 0.1 dex changes the oxygen deficit 
by 0.7 dex at a stellar mass of 10 9 M Q and by 0.15 dex at 
a stellar mass of lO 11 -/^©. This level of dust depletion in 
our model with a large oxygen deficit does not alter the 
conclusions of this paper. Given the large uncertainties 
on how dust depleted oxygen is in HII region and sys- 
tematic uncertainties of our metallicity calibration, we 
make no explicit correction for this effect. In Figure 12c, 
we vary the adopted metallicity by [-0.2, +0.1] dex. 

9.5. Abundance Gradients 

It is well established that star-forming galaxies in the 
local universe have spatial abundance gradients such that 
the central regions of galaxies are more metal rich than 
their outskirts (e.g. Zaritsky et al. 1994). In our esti- 
mate for the amount of oxygen in the gas-phase (see 
Section 8.2) we adopt an oxygen abundance determined 
from nebular emission within the 3 arc second SDSS fiber 
aperture. We determine the oxygen abundance of local 
galaxies from data selected to have a covering fraction of 
> 30% which is shown to be sufficient for reproducing the 
global abundance (Kewley et al. 2005). We reproduce the 
global luminosity weighted abundance whereas an unbi- 
ased method for estimating the gas-phase metal content 
of galaxies would be to use a global metallicity estimator 
that reflects the average total gas-phase abundance (i.e. 
a gas mass weighted abundance). However, given cur- 
rent observations it is not possible to assess what biases 
are introduced by adopting a luminosity weighted abun- 
dance for estimating gas-phase metal content. Surveys of 
abundance gradients in nearby galaxies currently under 
way along with a new generation of radio instruments 
capable of measuring the gas content of large samples of 
galaxies will likely resolve this issue. 

We speculate that the luminosity weighted oxy- 
gen abundance is likely to overestimate the gas mass 
weighted abundance. In star-forming galaxies where 
luminosity is dominated by young stars the luminos- 
ity weighted abundance reflects the abundance of star- 
forming gas. The Schmidt-Kcnnicutt relation for star 
formation relates the star formation surface density to 
the gas surface density by ^sfr oc Sg Q 4 s (Kennicutt 
1998b). Because of the non- linear relation between the 
two quantities, the luminosity weighted metallicity will 
not be consistent with the gas mass weighted metallicity 
but instead will be biased towards higher gas surface den- 
sity regions which tend to lie in the central, metal-rich 
regions of galaxies. Therefore, the average metallicity of 
the gas is likely to be lower than the inferred luminosity 
weighted abundance used in this study for estimating the 
oxygen mass in the gas-phase. However, given that obser- 
vations comparing the gas mass and luminosity weighted 
abundances are currently not feasible, we do not attempt 
to make a correction for this effect. 

10. DISCUSSION 

The census of oxygen for star-forming galaxies has been 
determined from the observationally constrained star for- 
mation and chemical histories. In Section 8 we derive the 



oxygen deficit determined by comparing the expected to- 
tal production of oxygen with what is currently found in 
the gas and stellar phases. We note that oxygen deficit 
represents the net deficit of oxygen accounted for by the 
stars and gas and excludes oxygen that may be expelled 
from the ISM and subsequently recycled back into the 
galaxies as some models predict (e.g Dave et al. 2011). 

One advantage to determining mass loss in the man- 
ner presented in this study is that a detailed physical 
mechanism governing mass loss is not required. The 
feedback processes thought to be responsible for mass 
loss have proven to be notoriously difficult to model and 
disentangling the various contributions observationally is 
challenging given that the physical mechanisms are not 
clearly understood (for review see Veilleux et al. 2005). 
Here we discuss some of the implications of our observa- 
tional constraints on mass loss and future prospectives 
for our models and results. 

10.1. Outflows 

Galactic outflows are a key physical process governing 
galaxy evolution (Veilleux et al. 2005). In the nearby 
universe they are commonly observed in starburst and 
post-starburst galaxies (e.g Rupke et al. 2005; Martin 
2006; Tremonti et al. 2007; Rich et al. 2010; Tripp et al. 
2011) and are ubiquitous in higher redshift star- forming 
galaxies (e.g Shapley et al. 2003; Weiner et al. 2009; Ru- 
bin et al. 2010; Stcidel et al. 2010). The outflowing gas 
and metals may remain gravitationally bound to galaxies 
in which case the material should be found to reside in 
the hot halos or it may escape the galaxy potential well 
altogether. 

The oxygen deficit may be partially resolved by oxygen 
present in the outer disks and hot halos of galaxies. Sev- 
eral studies have found flat abundance gradients in the 
outer disks of star- forming galaxies (Bresolin et al. 2009b; 
Werk et al. 2010, 2011; Bresolin et al. 2012). The amount 
of oxygen observed in the outer disks cannot be recon- 
ciled with the low levels of star formation and transport 
of metals from the inner disk is the most likely scenario 
explaining the flat abundance gradients though the phys- 
ical mechanism for the transport of metals to outer disks 
is not clearly understood. Oxygen has also been observed 
in absorption in the halos of star-forming galaxies (Chen 
& Mulchaey 2009; Prochaska et al. 2011). In a survey of 
42 galaxies conducted using the Cosmic Origins Spectro- 
graph onboard the Hubble Space Telescope, Tumlinson 
et al. (2011) report ubiquitous detection of OVI in the 
halos of star-forming galaxies. They conclude that the 
mass of heavy elements and gas in the circum-galactic 
medium (CGM) may far exceed whats found within the 
galaxies themselves. From the column density of OVI, 
they estimate that the CGM of 10 9 5 - 10 115 M Q stellar 
mass galaxies (out to 150 kpc) contains at least ~ 10 7 M Q 
of oxygen, where they have assumed an ionization correc- 
tion factor of 0.2 to account for oxygen in all ionization 
states. 

In Figure 14 we plot the oxygen deficit along with 
the lower limit mass estimate of Tumlinson et al. (2011) 
shown by the red dot-dashed line. At a stellar mass of 
10 9 ' 5 Mq, our model showing the largest oxygen deficit 
(solid curve) is consistent with the lower limit estimates 
of Tumlinson et al. (2011). These independent results 
rule out all except the highest oxygen deficit model (solid 
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Fig. 14. — The oxygen deficit with the lower limit of observed 
oxygen mass in the halos of star-forming galaxies from Tumlinson 
et al. (2011). The oxygen deficit is the same as in Figure 12d and 
the solid, dotted and dashed lines indicate the different oxygen 
deficit derived by varying the free parameters. 



line in Figure 14). A model with significant and ubiq- 
uitous outflows in star-forming galaxies is required for 
consistency with independent observations. As discussed 
in Section 7.4, our models suggest that in the case of 
significant outflows the oxygen deficit will scale (nearly) 
linearly with stellar mass. A linear scaling is consistent 
with the result found by Kirby et al. (2011) from their 
analysis of metal mass loss in Milky Way dwarf galaxies. 
At lO n M0, our model exceeds the lower limit estimate 
of Tumlinson et al. (2011) by ~ 2 orders of magnitude. 

Gas propelled to high velocities may escape from galax- 
ies all together. Estimates of the escape fraction of out- 
flowing material have been difficult to determine accu- 
rately mainly due to lack of constraints on halo drag 
(Veilleux et al. 2005). It may be that detailed estimates 
of the oxygen mass in the CGM of galaxies will account 
for the ~ 2 orders of magnitude greater oxygen deficit 
at M* ~ 1O U M implied by our models as compared to 
lower limit estimates of Tumlinson et al. (2011). How- 
ever, if the oxygen deficit is not resolved by oxygen found 
in the halos of galaxies, a natural reservoir for the re- 
mainder of the oxygen deficit would be the intergalactic 
medium (IGM) . This scenario is consistent with theoret- 
ical models where large scale galactic outflows transport 
metals from the galaxy to the IGM (e.g. Springel & Hern- 
quist 2003; Oppenheimer & Dave 2006; Kobayashi et al. 
2007; Oppenheimer et al. 2011). However, the kinemat- 
ics of oxygen in the CGM of local star-forming galaxies 
suggest that most of the oxygen is gravitationally bound 
(Tumlinson et al. 2011). 

A common method of estimating the escape fraction 
is to compare the outflow velocity to the escape ve- 
locity of the gravitational potential well of the galaxy. 
Early theoretical work suggested that mass loss occurs 
only in galaxies with shallow potential wells (e.g. Lar- 
son 1974; Dekel & Silk 1986; Mac Low & Ferrara 1999; 
Ferrara & Tolstoy 2000). However, recent observational 



and theoretical work has suggested that mass loss driven 
by AGN or star formation occurs in significantly higher 
mass galaxies as well (Strickland et al. 2004; Murray et al. 
2005; Weiner et al. 2009; Steidel et al. 2010; Fischer et al. 
2010; Feruglio et al. 2010; Sturm et al. 2011; Aalto et al. 
2012). 

It is not clear whether low-mass or high-mass galax- 
ies are responsible for the IGM enrichment. Though low 
mass galaxies may lose a larger fraction of their metals, 
they may not be a significant source of IGM enrichment 
due to their low rates of metal production (Martin et al. 
2002; Kirby et al. 2011). On the other hand, due to the 
deep potential wells of large galaxies, metals may not be 
efficiently ejected into the IGM. If the latter is true, our 
model would predict that the CGM and/or outer disks 
of massive star-forming galaxies (10 1:l Mq) should con- 
tain substantially greater quantities of oxygen than their 
lower mass counterparts (10 9 M Q ). Surveys of the CGM 
and studies of the outer disks of star-forming galaxies 
currently underway will provide important clues to re- 
solving the oxygen deficit implied by our models. 

10.2. Baryonic Mass Loss 

One possible interpretation of the oxygen deficit is that 
it has been driven out of the galaxy ISM by galactic 
outflows. If oxygen is expelled via outflows, we expect an 
even larger mass of gas expelled since outflowing gas will 
not be pure oxygen. We can estimate the total baryonic 
mass loss in two ways. First, we adopt an enriched wind 
model assuming that the metallicity of outflowing gas is 
equal to the Type II SNe oxygen yield (see Section 8.2). 
We adopt a value for P Q = 0.007 which is - 1.6Z Q . This 
value is consistent with estimates of Martin et al. (2002) 
for the metallicity of a starburst driven wind in NGC 
1569. Assuming outflows are enriched, we can estimate 
the total baryonic mass loss by the dividing the oxygen 
deficit by the oxygen yield. This is given by 



AM h 



M° 



M° - M° 



Po 



Po 



(29) 



AMb is the total amount of baryonic mass loss and AM 
is the oxygen deficit. 

We also derive an upper limit estimate of the total 
baryonic mass loss by assuming a uniform wind model. 
In the uniform wind model, the metallicity of the outflow- 
ing gas is the same as the ISM. In general, the metallicity 
of the outflowing gas is lower than the yield and therefore 
more gas is required to outflow in this model in order to 
yield the oxygen deficit derived in Section 7.4. In the uni- 
form wind model estimate P in Equation 29 is replaced 
by a SFR weighted metallicity, Z^. The SFR weighted 
metallicity and the total baryon loss are given by 
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and AMty > AM b . 

In the top panel of Figure 15 we plot the total baryonic 
mass loss as a function of stellar mass for galaxies in the 
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(Penton et al. 2004; Sembach et al. 2004), the majority 
of baryons are still missing. Some cosmological simula- 
tions favor the warm-hot intergalactic medium (WHIM) 
as the repository of the missing baryons (e.g Cen & Os- 
triker 1999; Dave et al. 2001; Cen & Ostriker 2006; Op- 
penheimer et al. 2011). However, observations of the hot 
gas at 10 5 — 10 7 K comprising the WHIM remain tenta- 
tive and no compelling evidence for the detection of this 
phase yet exists (Bregman 2007). The distribution of 
the WHIM material is unknown and hot halos of mas- 
sive galaxies are considered as a possible reservoir for 
substantial fraction of the missing baryons (Cen & Os- 
triker 2006; Tang et al. 2009; Kim et al. 2009), though 
this remains controversial (Anderson & Bregman 2010). 

An open question is what fraction of the missing 
baryons were ejected from galaxies through feedback pro- 
cesses and what fraction never accreted in the first place. 
We can compare our estimates of the total baryonic mass 
that can be associated with galaxies (baryonic mass plus 
the baryonic mass loss) with the expected baryon content 
of galaxies and their halos inferred from cosmological es- 
timates. The stellar-to-halo mass (SHM) relation param- 
eterizes the relationship between stellar mass of galaxies 
and the dark matter halos in which they reside. Moster 
et al. (2010) develop a statistical approach whereby ha- 
los and subhalos are populated in an N-hoAy simulations 
with the requirement that the observed stellar mass func- 
tion be reproduced. They parameterize the SHM as 



Fig. 15. — The top panel is the inferred baryonic mass loss and 
the bottom panel is the baryonic mass loss relative to the stellar 
and gas mass. The oxygen deficit is the same as in Figure 12d and 
the solid, dotted and dashed lines indicate the different oxygen 
deficit derived by varying the free parameters. The black and red 
curves are determined by adopting an enriched outflow (Equation 
29) and uniform wind model (Equation 31), respectively. 

local universe. The black and red curves are estimates 
adopting an enriched outflow (Equation 29) and uniform 
wind model (Equation 31), respectively. In the bottom 
panel of Figure 15 we plot the logarithm of S b which we 
define as 

AM b 

5b= M^Ml- (32) 

5b is the ratio of the total baryonic mass loss to the bary- 
onic mass of the galaxy, where the baryonic mass of the 
galaxy (gas + stellar). St represents the ratio of the mass 
of gas to be cycled in and out of the galaxy compared to 
the current baryonic mass of the galaxy. 

In the currently accepted cosmological model ( ACDM) , 
the universal baryon fraction is precisely determined 
from the cosmic microwave background, the observed 
baryon acoustic oscillations and the Hubble constant. 
The universal baryon and dark matter density revealed 
by these observations is given by fit, = 0.0456 ± 0.0016 
and f2 c = 0.227 ± 0.014, respectively (Komatsu et al. 
2011). Studies attempting to account for the baryons 
find that only a fraction of the expected baryons are ob- 
served in the low-redshift universe (Fukugita et al. 1998; 
Fukugita & Peebles 2004; Nicastro et al. 2005; Sommer- 
Larsen 2006; Shull et al. 2011). 

Only about a tenth of the baryons are found in the stars 
and gas of galaxies (Bell et al. 2003a). While the Lya 
forest at low redshifts can account for another ~ 30% 
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Here Mh is the halo mass and (j^J , Mi, j3 and 7 are 

free parameters. The relation evolves with redshift and 
the parameters are given by 
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We can estimate the expected baryon content of galaxies 
from the universal ratio of baryonic to dark matter given 

by f bc = n b /n c ===== 0.201 ± 0.014. 

Figure 16 demonstrates the missing baryon problem 
for galaxies. The dashed cyan line is the current total 
baryonic (gas + stellar) mass content of galaxies plotted 
as a function of stellar mass. The dot-dashed blue line 
is the expected total baryonic mass content of galaxies 
assuming the universal baryonic to dark matter ratio and 
the SHM relation of Moster et al. (2010) and is given by 
M b = f bc Mh ■ The solid black and red curves are the sum 
of the baryonic mass and the baryonic mass loss using the 
two estimates given by Equations 29 and 31. The solid 
black and red curves can be interpreted as estimates of 
the baryon content associated with galaxies and are the 
amount of baryons currently found in local star-forming 
galaxies plus what was once in the galaxies but has since 
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The quantity AM^ is the total amount of baryonic mass 
loss and M»/(l — i?) represents the total amount of star 
formation. In the literature the instantaneous mass load- 
ing factor is defined as ry, = M w /M*. The instantaneous 
mass loading factor is the mass loss rate divided by the 
SFR. The effective mass loading as we have defined it in 
Equation 35 is the instantaneous mass loading factor, rji, 
averaged over the star formation history. In Figure 17 
we plot the effective mass loading factor as a function of 
stellar mass. Here, we have adopted R — 0.5. Assuming 
a smaller value for R would lower the estimate. Even 
adopting our upper limit estimates of baryon mass loss, 
the effective mass loading for star-forming galaxies in the 
local universe is < 1. 



Fig. 16. — The dashed cyan curve is the total baryonic (gas + 
stellar) mass of galaxies in the local universe. The dot-dashed 
blue line is the inferred baryon content from cosmological fraction 
and is given by /j, c . The solid black and red lines are the total 
baryonic mass plus the baryon mass loss for an enriched wind model 
(Equation 29) and uniform wind model (Equation 31), respectively. 
In our estimate we have adopted the oxygen deficit given by the 
solid curve in Figure 12d. 
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Fig. 17. — The effective mass loading parameter given in Equa- 
tion 35 plotted against stellar mass. The red and black curves 
give the effective mass loading factor for our upper and lower limit 
estimates of the total baryon mass loss, respectively. 

been cycled out through outflows. 

The baryon content associated with galaxies is still 
substantially lower than the inferred content from cos- 
mology. The straightforward interpretation of Figure 16 
is that the missing baryons were never accreted onto lo- 
cal star-forming galaxies and unless a large reservoir of 
gas is found in the halos of star-forming galaxies it is 
likely that the missing baryons reside in the IGM. Using 
observationally motivated constraints for the mass and 
radii of hot halos, Anderson & Bregman (2010) come to 
a similar conclusion. 

From our estimate of the baryon mass loss we can de- 
rive an effective mass loading factor, 77, which is given 



10.3. Future Prospects 

A complete and self-consistent theory of galaxy evolu- 
tion will require a detailed account of galaxy growth and 
chemical evolution along with physical mechanisms gov- 
erning these processes. In this contribution we present 
empirical models attempting to self-consistently inte- 
grate chemical evolution and galaxy growth. Our self- 
consistent approach is complementary to cosmological 
simulations and semi-analytical models and the self- 
consistent census approach presented here should be 
compared with those models. Both theoretical and obser- 
vational advances are crucial to constraining the model 
results developed in this study. The method of analysis 
used in this study provides useful tests for consistency 
of a diverse set of observations and theories related to 
chemical properties of galaxies. Here we address future 
prospects for improvement. 

One of the greatest outstanding astrophysical prob- 
lems is the large uncertainty in the absolute calibration 
of the nebular abundances scale. High S/N observations 
of large sample of HII regions covering a broad range of 
physical parameters will be extremely important in sta- 
tistically establishing and testing diagnostics. Such ob- 
servations will be crucial in developing empirical calibra- 
tions relying on recombination lines which are thought 
to be less susceptible to effects of temperature variations 
associated with auroral lines (Esteban et al. 2002; Peim- 
bert & Peimbert 2005; Bresolin 2007). A complementary 
approach will be to use the recently developed wide-field 
integral field spectrographs to observe nearby HII regions 
in order to understand the discrepancy between empiri- 
cally and theoretically calibrated strong-line methods. A 
well calibrated diagnostic applied to large data sets inves- 
tigating the MZ relation at higher redshifts and to lower 
stellar masses will provide important constraints for our 
understanding of chemical evolution and our census of 
oxygen in the stellar and gas-phase. 

Theoretical models incorporating stellar rotation and 
mass loss into computations of nucleosynthesis for stars 
at all masses are yet to be developed (Romano et al. 
2010). Current models of stellar nucleosynthesis are not 
yet able to explain the full diversity of chemical abun- 
dance patterns observed. While the oxygen yields are 
much better constrained by models owing to its primary 
origin, there is still a factor of ~ 2 discrepancy. A sin- 
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gle, self-consistent model of nucleosynthesis able to re- 
produce the abundance ratios observed in the Milky Way 
and other nearby galaxies will likely alleviate some of the 
discrepancy, thus constraining the total oxygen produc- 
tion in star-forming galaxies. 

Part of the uncertainty in chemical evolution models 
rests on the assumption of a particular IMF (Romano 
ct al. 2005). Adopting a constant IMF is problematic for 
galaxies, where even a constant stellar IMF, could lead 
to variations in the integrated galactic IMF (Kroupa & 
Weidncr 2003) and studies of the IMF from integrated 
measurements of star-forming galaxies indirectly indicate 
variations (Hovcrsten & Glazebrook 2008; Lee et al. 2009; 
Meurer et al. 2009). A varying IMF would have im- 
portant implications for chemical evolution models (e.g. 
Romano et al. 2005; Koppen et al. 2007; Calura et al. 
2010). Despite the mounting evidence, no direct evidence 
of IMF variations in star-forming galaxies is currently 
available leading to a lack of consensus on the constancy 
of the IMF. Distinguishing whether a universal IMF or 
IGIMF formulation provides a better description of large 
scale star formation will be important for chemical evo- 
lution studies of galaxies. Large statistical studies of 
the integrated properties of galaxies observed over cos- 
mic time, in particular observations of the FUV and FIR 
properties which are now becoming available, will be use- 
ful in establishing any systematic variations of the IMF. 
In our study, tighter constraints on the IMF will alleviate 
much of the uncertainty in the total amount of oxygen 
produced and the return fraction. 

Total gas masses are required to measure the abso- 
lute content of metals within galaxies from observed rel- 
ative abundances. In the local universe, large surveys of 
atomic and molecular gas have been conducted (Heifer 
et al. 2003; Walter et al. 2008; Lcroy et al. 2009; Sain- 
tonge et al. 2011). Molecular gas at higher redshifts 
has also been detected in star-forming galaxies (Daddi 
ct al. 2008, 2010; Tacconi et al. 2010). The next gen- 
eration of radio and sub-mm observatories such as the 
Atacama Large Millimeter/submillimeter Array and the 
Square Kilometer Array will revolutionize the study of 
cold gas in the universe allowing us to probe larger sam- 
ples to far greater depths. Currently, we are only able to 
estimate the total oxygen in the gas-phase of local star- 
forming galaxies. Measurements of cold gas in nearby 
and distant star-forming galaxies taken together with a 
well calibrated gas-phase metallicity diagnostic will allow 
us to track the total mass, not just relative abundance, 
of oxygen in the gas-phase of star-forming galaxies over 
cosmic time. 

The structure, content and chemical composition of 
outflowing gas is highly uncertain. Understanding the 
physical properties of galactic winds is crucial for es- 
timating the total mass of gas outflowing from galax- 
ies. Observations reveal that galaxy scale outflows have 
a complicated multiphase structure (see Veilleux et al. 
2005) and observations suggest that structure and com- 
position vary in each of the phases (e.g. Tripp et al. 2011). 
Multi-wavelength observations are required. The cold 
gas component can be observed in absorption lines of 
neutral and low ionization state metals and the current 
generation of integral field spectrographs on 8-10m class 
telescopes and wide-field integral field spectrographs on 
4m class telescopes will provide census of galactic scale 



outflows in local galaxies and presence of galactic winds 
in the distant universe. Hotter phases can be observed 
using space-based UV spectrometers such as the Cosmic 
Origins Spectrograph (COS) onboard Hubble. A crucial 
but inaccessible phase is the so-called wind fluid which 
drives the stellar winds. A hard x-ray telescope with 
high spatial resolution and sensitivity is required to ob- 
serve this phase. For the foreseeable future, astronomers 
will likely have to rely on detailed theoretical models and 
indirect observations to understand the wind fluid. 

COS will also provide important insight into the 
circum-galactic and warm-hot ionized mediums, both 
thought to be important repositories of baryons. Studies 
of these regions will begin to reveal the baryonic con- 
tent of the hot haloes and intergalactic medium sur- 
rounding galaxies. Observations of these regions along 
with a census of baryons associated with galaxies will 
be crucial to potentially identifying the large fraction of 
missing baryons in the local universe and resolving this 
long-standing problem (Brcgman 2007). A benchmark 
for these studies is to resolve the oxygen deficit in star- 
forming galaxies presented here. 

In this study we have applied the best theoretical and 
observational constraints available in undertaking a cen- 
sus of oxygen in star-forming galaxies. We are unable 
to draw any strong quantitative conclusions from the 
models developed owing to the large uncertainties as- 
sociated with both our adopted model parameters and 
observational inputs. Nonetheless, this study represents 
one of the first attempts to self-consistently account for 
the stellar mass growth and chemical evolution of galax- 
ies. As theoretical and observational advances allow for 
ever greater constraints, we hope the approach taken in 
this study will prove to be an important ingredient in 
testing and developing a fully self-consistent theory of 
galaxy evolution. 

11. SUMMARY 

In this study, we have consistently incorporated chem- 
ical evolution in the framework of stellar mass growth in 
order to conduct a census of oxygen in local star-forming 
galaxies. We are able to estimate the total oxygen pro- 
duction from total amount of star formation inferred 
from the current stellar mass. The mass of oxygen in 
the gas-phase is constrained by the observed MZ rela- 
tion in the local universe and the relation between gas 
fraction and stellar mass. The most difficult to constrain 
observationally is the amount of oxygen locked up stars. 
Our empirical models self-consistently incorporate stel- 
lar mass growth and chemical evolution, thus allowing us 
to track the metallicity of the gas from which stars are 
formed over cosmological timescales and giving us empir- 
ical constraints on the oxygen mass locked up in stars. 
The main results of this study are given below. 

1. We conduct a census of oxygen and show that the 
amount of oxygen in the stellar and gas phase 
of galaxies does not fully account for the total 
amount of oxygen produced. We conclude that the 
most straight-forward interpretation of the oxygen 
deficit is that oxygen has been expelled from galax- 
ies by outflows. Our results establish the need for 
ejective feedback in normal star-forming galaxies. 
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2. We compare our oxygen deficit with the observed 
lower limit of oxygen found in the CGM of star- 
forming galaxies and conclude that oxygen mass 
loss is a ubiquitous process in star-forming galax- 
ies. Furthermore, the oxygen deficit in our pre- 
ferred model scales with stellar mass and we predict 
that cither more massive galaxies should be found 
to contain a greater mass of oxygen in their halos 
or that oxygen escapes the galaxy potential well al- 
together and therefore massive galaxies contribute 
for IGM enrichement. 

3. We estimate the total amount of mass lost from the 
ISM of star-forming galaxies and find that it is a 
small fraction of the total baryon content expected 
from the cosmological baryon density. We conclude 
that only a small fraction of the total baryons in the 
universe ever cycled through star-forming galaxies. 

Our empirical model provides an important test of self- 
consistency for many physical processes governing galaxy 
evolution. Future theoretical and observational advances 
will provide ever increasing constraints on the census of 
oxygen in star-forming galaxies and our models provide 
important benchmarks with which to compare theory 
and observation. 
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APPENDIX 
Al. Time Dependent Mass Return 

Here we test the validity of our assumption of a con- 
stant return rate with instantaneous recycling. Jung- 
wiert et al. (2001) develop an algorithm to account for 
continuous, time dependent mass loss. They give the gas 
return rate as 



M R (t) = / ¥(*')/„ 
Jo 



.(t-t')dt'. 



(36) 



Here Mfl(t) is the gas return rate as a function of time 

and f m r(t — t') is the mass loss rate at time t for a single 
burst stellar population. The mass loss then is given as a 
convolution of the SFR with the mass loss rate (Equation 
36). The gas return for various IMFs is well parameter- 



ized by the equation 

/ ror (t) = <7 ln(y + lY (37) 
Here both Co and A are constants depending on the par- 
ticular choice of IMF (values are given in Table 1 of Leit- 
ner & Kravtsov 2011). We determine the gas return rate 
using a Chabrier and a Chabrier steep IMF. These two 
IMFs have been chosen because they are representative 
of the upper and lower mass loss rates of most IMFs (see 
Figure 1 of Leitner & Kravtsov 2011). 

We employ the similar algorithm for determining a self- 
consistent gas return rate and stellar mass history as out- 
lined in Leitner & Kravtsov (2011). We determine a star 
formation history assuming Equation 19 with the return 
fraction R — 0.35. Using this star formation history we 
then determine the gas return rate from Equation 37. We 
then replace the constant return fraction with the time 
dependent return rate in Equation 18 such that we have 



dM* 
dz 



R- 



(38) 



We integrate this equation to determine a new star 
formation history which we then use to redetermine a 
new gas return rate via Equation 36. We iterate un- 
til AM* < 0.01, where AM* is the change in the final 
stellar mass between iterations. This procedure typically 
takes < 10 iterations. In this procedure, the gas returned 
has the metallicity of the natal gas, thus relaxing the as- 
sumption of instantaneous recycling. 

Following the procedure outlined in Section 8, we de- 
termine the oxygen deficit using the time dependent re- 
turn rate formulation. The oxygen deficit is shown in 
Figure 18 where we have used P Q — 0.008 and AZ g = 0. 
The black curves are the same as Figure 12d and are 
determined using a constant, instantaneous return rate. 
The solid and dashed red curves are the oxygen deficit 
determined from the Chabrier and Chabrier steep IMFs, 
respectively. This figure convincingly demonstrates that 
a constant return rate with instantaneous recycling is a 
valid simpl ification in . qur census of oxygen , 
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Fig. 18. — The black curves are the same as in Figure 12d. The 
red curves are the oxygen deficit determined by applying a time de- 
pendent return rate using the Chabrier (solid curve) and Chabrier 
steep (dashed curve) IMFs. 
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